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ABSTRACT 

Thirty- one high dispersion coude spectrograms of Q Oph 
(V = 2. 57, MK-type 09. 5V, b H = +23°) and seven of Q Per 
(V = 2. 83, MK-type Bllb, b = -17°) have been numerically syn- 
thesized to produce high resolution, low-noise spectra in the in- 
terval AX 3650 - 4350 A that yield more accurate data on atomic and 
molecular absorption in well-defined regions of the interstellar medium. 
The detection threshold is improved by as much as a factor 5 
over single plates. Several new interstellar lines have been 

discovered in the Q Oph - 15 km sec ^ cloud tand the Q Per + 13 km 

-1 13 + « , o 

sec cloud. C H A4232 ft. and KI AA4043 and 4047 A are seen in 

Q Oph for the first time anywhere in the interstellar medium, and 

the previously observed Ca I A4227 A and Fe I XA3720 and 3860 A 

for the first time in both Q Oph and Q Per. Their equivalent widths 

range from 0. 3 to 2 mA. Ca I X.4227 A, with the strong Ca II H and K 

interstellar lines, yields the ionization equilibrium and hence the 

electron density. Reasonable limits on the photoionization rate 

computed from recently proposed galactic and stellar radiation 

. 3 _ 1 

fields yield n = 0.16 — 0.55 cm in the c. Oph -15 km sec 
e 

- 3 - 1 

cloud and n = 0. 044 — 0. 14 cm in the C Per +13 km sec cloud. 

e ■■ 

‘ I ; r ' Z 
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-3 

The neutral hydrogen density is estimated to be 50 — 7 0 0 cm in 

_3 

the £ Oph cloud and 10 — 80 cm is £ Per for reasonable assump- 
tions of the cosmic ray ionization rate of hydrogen and the abundance 
of easily ionized trace elements. Sodium and potassium appear in 
the Q Oph cloud approximately at their solar system abundances 
relative to hydrogen, but calcium, iron, and aluminum seem to be 
deficient by large factors — 2500, 300, and > 100, respectively. 

The abundances in £ Per relative to sodium are similar, suggest- 
ing that most of the iron, aluminum, and perhaps other metals is 
bound up in the interstellar grains. 

C"^H + \4232 A in £ Oph yields a [ ratio of 81 

(+ 113, -16), consistent with terrestial value 89. This suggests 
13 

that C -rich stars have not contributed much to the interstellar 

medium, at least in the vicinity of the sun, and that its compo- 

9 

sition has not changed much in the past 5 X 10 years. 

The CN, CH, and CH + bands at W3875, 4300, and 4232 A 
provide in all one precise measurement of the microwave background 
radiation brightness temperature and three upper limits at milli- 
meter and submillimeter wavelengths. The CN excitation tempera- 
tures, corrected for local processes, yield T (\2. 64 mm) = 

2. 75 ± 0. 10° K in Q Oph and 2. 66 ± 0. 12° K in £ Per. The upper 


iv 


limits on the background radiation intensity in units of erg cm 

-1 -1 -1 -14 -14 

sec ster Hz are 0. 72 X 10 at 1. 32 mm (CN), 1. 66 X 10 

at 0. 56 mm (CH) and 3. 75 X 10 " 14 at 0. 36 mm (CH + ). The CN 
excitation in six stars, including late B and early A giants, is 
found to be consistent with 2. 8° K and consistent with a universal 
CN excitation mechanism. 
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I. INTRODUCTION 

Soon after the discovery of the microwave background 
radiation by Penzias and Wilson (1965), -it was recognized that 
the optical absorption lines of the interstellar molecules CN, CH, 
and CH + might be used to determine the radiation brightness 
temperature at certain millimeter and submillimeter wavelengths 
(Field and Hitchcock 1966; Thaddeus and Clauser 1966). The back- 
ground photons at the wavelengths of molecular resonances will 
rotationally excite the molecules in interstellar space, and hence 
will be detectable optically by the presence of band spectra. The 
relative intensities of absorption lines which arise from different 
lower levels will yield the brightness temperature of the background 
radiation at the transition wavelength connecting those levels. 

It had long been known that the first excited rotational 
level of CN, whose J = 0 -* 1 transition resonance falls at 
\ 2. 6mm, was populated in interstellar space. On the basis 
of visual estimates by Adams (1941) of the intensities of the 
R(0) and R(l) lines of the interstellar (0,0) violet band in the 
bright O star Q Oph, McKellar (1941) computed the excitation 
temperature of this level to be 2. 3° K. Hardly any significance 
was attached to this result at the time, it presumably being 
thought that the excitation was due to "local" interstellar 
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processes. 
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The first modern use of CN as a thermometer at the 
\ 2. 6mm resonance yielded excitation temperatures against Q Oph 
and £ Per that agreed both with the Penzias and Wilson (1965) 

\ 7. 35cm radio measurement and with each other (Field and 
Hitchcock 1966; Thaddeus and Clauser 1966). The agreement 
between the two molecular excitation temperatures suggested 
that the effect of local interstellar excitation mechanisms was 
small and that the excitation temperatures were actually quite 
close to the brightness temperature. 

Interstellar lines observed optically are never strong, and 
those of CN in particular are found only in a handful of stars. 

The technique of plate synthesis, where a number of spectro- 
grams of a star are added together, may be used to improve 
the signal- to- noise ratio of the weak molecular bands. The 
first application of this technique to background radiation mea- 
surements made use of existing plates, many obtained by Adams 
and Dunham in the 1930's and 40' s (Clauser 1970; Clauser and 
Thaddeus 1972). In addition to more accurate \ 2. 6mm CN 
temperatures in Q Oph and £ Per, this work established upper 
limits at three shorter wavelengths — \\ 1. 32, 0. 56, and 0. 36 
mm of CN, CH, and CH + , respectively. Data from a number 
of other stars were fcited as evidence for the directional 
invariance of the CN excitation. 
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This dissertation began as an attempt to still further 
improve the observational data on CN, CH, and CH + by the 
acquisition of new observations. About fifty high resolution 
plates were obtained at the Lick Observatory 120- inch telescope 
coude spectrograph, including 31 of the most favorable star, 
q Oph, and 7 of Q Per, where the CN lines (but not those 
of CH or CH + ) are about as strong as in Q Oph. Preliminary 
results of the Q Oph plate synthesis have been reported 
(Bortolot, Clauser, and Thaddeus 1969). In this dissertation, 
a significant improvement over those findings will be presented, 
along with new data on Q Per and other stars which pertain 
to the question of directional invariance. 

A good part of this dissertation will treat the syntheses, 
covering the entire interval XX3650 to 4350 A, of Q Oph and 
£ Per, both stars exceptionally rich in interstellar molecular 
and atomic absorption in spite of their relative closeness to 
the sun and moderate reddening. Although previous work on 
the synthesis technique involved only narrow intervals about 
strong absorption features (e. g. Augason and Herbig 1967; 
Bortolot et al . 1969), much effort was expended in the present 
work to synthesize entire plates, to take advantage of the con- 
siderable information on the interstellar medium they contain. 

10 



The great increase in sensitivity of the syntheses over single 
plates (more than five-fold for Q Oph) permits more accurate 
data to be obtained on interstellar atomic and molecular abun- 
dances than previously possible. These abundances in turn will 
form the basis of more accurate calculations of the ionization 
equilibrium and molecular processes. 

The interstellar line profiles in the spectra of the two 

f* 

stars we consider here are reasonably simple — that is, 
resolvable into distinct clouds with little overlap. We thus 
avoid the problem of interpretation occurring in studies of dis- 
tant stars where the lines are formed in multiple overlapping 
clouds. The analysis of such lines must lead to some sort of 
average physical conditions which conceivably might not apply 
to the average cloud. 

In the next section we will deal with the Q Oph observations, 
first obtaining the atomic and molecular column densities, then 
deriving from them the ionization equilibrium, and finally deducing 
the abundances. Of considerable interest will be the upper limits 
on undetected lines. The third section, on Q Per, will present 
essentially this same analysis on a different cloud. The mea- 
surements of the background radiation and a discussion of the 
molecular excitation mechanisms in the Q Oph cloud will comprise 
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II. THE -15 KM SEC" 1 CLOUD IN Q OPHIUCHI 

The interstellar line spectrum of Q Oph (V = 2. 57, MK-type 
09. 5V, b 11 = + 23°) has probably been more thoroughly studied 

than that of any other star. Adams (1949) noted there the presence 
of every then-known interstellar line with the exception of those 
of Fe I and Ca I. The optical lines in Q Oph appear to be resolvable 
into distinct clouds which hardly overlap. At conventional dis- 
persions, they are double, the strongest component lying at 
-15 km sec 1 and a weaker one at -29 km sec 1 . The high 
resolution Fabry-Perot observations of Hobbs(1965) have in- 
dicated structure at the bottom of the strong -15 km sec 1 line 
and at least four other components, near -9, -26, -28, and 
-33 km sec 

Herbig (1968) has given a quite thorough discussion of the 
physical conditions and the atomic and molecular abundances 
in the -15 km sec 1 cloud, and it would be worthwhile to briefly 
summarize his results. A curve of growth analysis has yielded 
doppler constants for Na I, Ca II, and CH + , and the column 
densities of Na I, Ca II, K I, CN, and CH. Although there is 
an extensive H II region surrounding Q Oph, roughly 15 pc 
in radius, it is argued that the -15 km sec 1 cloud is wholly an 
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H I region on the grounds that (1) the pattern of the optical 
absorption lines is consistent over an area of the sky considerably 
larger than that covered by the ionized region; (2) the 21-cm 
line velocities show general agreement with this pattern; (3) if 
the Na I lines were formed in the H II region, the Na/H abundance 
ratio would.be inconsistent with the solar system value; and 

(4) the rms turbulent velocities observed in H II regions are an 

0 

order of magnitude greater than those determined from the curve 
of growth analysis in the -15 km sec cloud. As the H II region 
is expected to be at least roughly sperical in shape, Herbig placed 
the cloud at a minimum distance of 15 pc from Q Oph. On the 
assumption that Q Oph dominates the radiation field and that 
there is no attenuation by the cloud, he computed the electron 
density in ionization equilibrium (see Equation (1), below) at a 

point 20 pc from the star. From an upper limit on the Ca I/Ca II 

-3 

line ratio, he obtained n < 1. 9 cm , but on the assumption that 

the Na/H abundance ratio had the solar system value 2. 4 X 10 ^ , 

• —3 

he deduced the ratio Na I/Na II and found n = 0. 35 cm . With 

e 

this electron density and the assumed radiation field, potassium 
was found to be normally abundant relative to sodium, while 
calcium was deficient by more than three orders of magnitude 
and titanium by at least two. He also deduced that the cloud was 

i JL3 

6 
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i 
! 
j 

j 

possibly a thin (0. 15 pc) sheet with the rather high neutral 

-3 

hydrogen density 500 — 900 cm 

In this chapter we will apply the plate synthesis technique 
to improve Herbig's observational data in the spectral region 

o 

M3650 — 4350 A, and to search for new interstellar lines. 

We then will discuss the physical conditions and abundances 
in the cloud, taking advantage of recent observational and 
theoretical work not available to Herbig, particularly in the 
discussion of the radiation field and the neutral hydrogen 
density. 

a) Spectrograms 

A total of 31 spectrograms of Q Oph in the region AA3650 — 

43 50 A were obtained with the coude spectrograph of the Lick 

Observatory 120-inch telescope in 1968-9. The dispersion with 

= -l 

the 160-inch camera was 1. 3 A mm and the resolving power 
A/AA was 85 000. The spectrograms were 20 inches long, 
consisting of two 10 -inch plates exposed end-to-end in the spect- 
rograph plate holder. The emulsion was Eastman Kodak Ila-O, 
baked at 50° C for 60 hours, and developed in D-76 for 15 min- 
utes. The projected slit width was 35 microns and the spectra 

14 
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were quite wide — about 5 mm, in order to store a good deal of 
information on each plate. The average exposure time was 
45 minutes. 

Densitometry was carried out on the D. W. Mann model 
1140 microdensitometer at the Goddard Institute for Space Studies. 
The output from the photomultiplier of the densitometer was 
digitized and, together with the comparator screw position, was 
recorded on magnetic tape every 4 microns of plate travel. 
Subsequent processing of the spectrograms was entirely digital. 

The methods of wavelength and intensity calibration and digital 
spectrum synthesis are described in Appendix A. The equivalent A 
widths were measured numerically: an intensity continuum was 

first fit to the vicinity of the line; the spectrum was then normalized 
to this continuum; finally, a numerical integration over the line 
gave the equivalent width. The procedure for numerical deter- 
mination of upper limits is discussed in Appendix B. 

b) Results of Plate Analysis 

An example of the marked increase in sensitivity provided 
by the plate synthesis of Q Oph is shown in Figure 1. The Fe I 
line at X3860 A shown in this figure cannot be seen on a single 



good plate, is barely detectable in a synthesis of six plates, but 

is unmistakable when all 31 plates are synthesized. The equivalent 

0 

width of this line is only 1 mA, or about one-half that of the faintest 
lines discernable on a single high quality plate. 

The atomic and molecular lines detected in the range 
A.A.3650 — 4350 A are shown in Figures 2 and 3. Table 1 lists the 
measured equivalent widths and radial velocities, and upper limits 
to the equivalent widths of the most important lines not detected. 

Herbig's (1968) data are also included for the sake of comparison. 

/ 

1 o | 

Three lines, K I \4044 A, X4047 A, and C H A4232 A, are 

detected in the interstellar medium for the first time. Three other 

Ca I A.4227 A and Fe I A. 3720 A and \ 3860 A are new in Q Oph. 

These lines, as Figures 2 and 3 show, fall at -15 km sec \ the 

velocity of the strong atomic and molecular lines. The Ca I line, 

which with H or K of Ca II allows determination of the ionization 

2 

equilibrium, has been seen previously against y Ori and several 

13 + 

other stars. A preliminary analysis of the C H line has been 
given (Bortolot and Thaddeus 1969). 

A number of interstellar lines were sought but not found. 

A search for all the resonance lines of the elements tabulated in 
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the Charlotte Moore Tables (1945) falling in the interval XX 3650 — 

4350 A yielded upper limits of 0. 5 mA or less, with the upper limit 
corresponding to a confidence level of 99 percent. Accordingly, 

o 

0. 5 mA may be taken as a general upper limit for any missing line. 
Some upper limits of particular importance are included in Table 1, 
and will be. discussed in later sections. 

In order to determine the optical depths and column densities 
of the atoms and molecules observed, some assumption concerning the 
curve of growth is necessary. . J3y fitting the equivalent widths of ■ 
the Na I and CH + lines to a curve of growth for a gaussian velocity 
distribution, Herbig. (1968) obtained the doppler constants b(Na I) = 

i 

2.43 km sec and b(CH + ) = 0. 85 km sec . Uncertainties in the 
M + 

CH Franck-Condon factors make the latter value somewhat 
doubtful. Hegyi, Traub, and Carleton (1972, to be published) have 
recently resolved the CN R(0) line in Q Oph with a Fabry-Perot 
spectrometer. When their linewidth is corrected for instrumental 
broadening and line saturation, b(CN) is found to be 1. 36 km sec \ 
This doppler constant, if caused by thermal line broadening, 
would indicate a kinetic temperature in excess of 1000° K— far 
above that of normal H I regions — and hence the line broadening 
is probably caused by microturbulence. On the assumption that 
the CN, CH, and CH + molecules which are observed are well mixed 

17 
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together, it will be assumed that b = 1. 36 km sec applies to all 
the molecular lines. Fabry-Perot work in progress by Hobbs (1972, 

-j- o 

private communication) on CH \4232 A will provide a check of 
this value. The doublet ratio method of Stromgren (1948) is well 
suited for determination of b(Ca II) from the moderately strong H 
and K lines and yields b(Ca II) = 1. 92 km sec \ The Ca II column 
density and those of the molecules, determined with the doppler 
constants just discussed, are given in Tables 2 and 3: the column 
densities in these tables corresponding to the weak, presumably 
unsaturated, lines are taken from the linear part of the curve of 
growth which is independent of doppler constant. 

The CN, CH, and CH + lines near U3875, 4300, and 4232 A 
used for determination of the microwave background temperature 
are shown in Figures 4, 5, and 6; Table 4 contains the equivalent 
widths and corresponding optical depths of these lines. Upper 
limits on the equivalent widths of those lines not observed were 
calculated according to the procedure of Appendix B. For these 
lines, the intensity continuum was taken to be the best straight line 

o 

fit to a 1 to 2 A interval about the expected wavelength; this is shown 
superimposed on the lowest spectrum in Figures 4 through 6. 

It is well known that systematic differences of the order of 

25 percent exist in equivalent width measurements at different 

18 
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observatories (see, e. g. , Wright 1969, Frisch 1972), but for 
determination of excitation conditions and the microwave back- 
ground temperature only the ratio of line intensities is important. 

With respect to relative intensities, it is estimated that the errors 
in intensity calibration (see Appendix A) are negligible in comparison 
to the statistical uncertainties of equivalent width measurement, at 
least for faint lines. Statistical uncertainties are discussed in 
Appendix C, and are given for the CN lines in Table 12. 

’ ' ' ... p . . 

Table 1 shows that no systematic difference in radial velocity 
exists between various atomic and molecular lines. This obviously 

i 

supports the idea that the -15 km sec ^ cloud is rather homogeneous, j 

i 

j 

and that the atoms and molecules are all well mixed. A possible 

-j- o o 

exception to this is found in CH X3958 A, which is about 0. 01 A to 

the blue of its expected wavelength. Herbig seems to observe the j 

same effect, and it appears also in the Q Per observations to be 
discussed. The most likely explanation is that the rest wavelength 

i 

of Douglas and Herzberg (1942), which was not directly measur ed ; 

but was computed from the molecular constants derived from other 
bands, is somewhat in error. 

i 

c) Electron Density and Atomic Abundances 

We turn now to the problem of extracting from the column 
densities as much information as possible on the physical conditions 
and abundances within the -.15 km sec ^ cloud. 

19 
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i) Ionization Equilibrium 


-1 


The detection of Ca I X4ZZ7 A in the -15 km sec cloud permits 
there for the first time a direct evaluation of the ionization equilibrium 
and the electron density. In ionization equilibrium (see, e. g. , Spitzer 
1968), 


n ( X r +1 ) n e b 

n(X r ) ' 


(1) 


where n(X^) is the number density of the element X 
th 

in the r stage of ionization, n is the electron density, F is the 

e r 

photoionization rate, and is the total recombination rate to the 
th 

r stage of ionization. The photoionization rate is given by 


r = h 1 J a u \d\ 


X X 


sec 


-1 


( 2 ' 


where a is the photoionization cross section, u is the radiation 

A. A. 

density, and h is Planck's constant (we have dropped the subscript r 
in Equation (2) since we are dealing only with neutral and singly 
ionized atoms). Thus, the observed Ca 11/ Ca I density ratio per- 
mits computation of n^ as a function of the radiation field; in turn, the 
electron density allows one to compute abundances of the elements 
observed, provided that their photoionization cross sections and 
recombination rates are known for all accessible stages of ionization. 
The only dependence of Equation (1) on the kinetic temperature is 


m 


a 



Th i nly d c p c rt-d c n -Q e of E g jiiilJ-aa- ( ' 4 -en t Ttg~icrr 


7-rf-K^-te«^e^3±nT^' is 

-0 4 

through a ~ T . Estimates of T ^ range from 50 - 125° K and 
100° K is adopted. The resulting uncertainty in the ionization equi- 
librium is only about 25 percent - at most comparable to that of 
other quantities in the equation - and will be neglected. 

ii) The Radiation Field 

The radiation field of optical photons in the -15 km sec ^ cloud 
is presumably the general galactic field plus that of Q Oph, attenuated 
by the cloud itself. Instead of the galactic spectrum computed by 
Zimmermann (1964) which Herbig used in his analysis, we will 
adopt the field recently calculated by Habing (1968) on the basis of 
new measurements of the interstellar extinction and the ultraviolet - 
flux from OB stars. Using rocket, satellite, and ground-based 
data to separate the airglow and zodiacal light contributions to the 
sky brightness, Lillie (1968) has found a radiation density of 

-18 -3 o -1 

35 — 55 X 10 erg cm A at A.2000 A, in good agreement with 
Habing's values 30 X 10 erg cm ^ ^ at \2200 A and 50 X 10^ 

-3c -l 

erg cm A at A. 1400 A. For the stellar field we will adopt, as 
did Herbig, the radiative flux from the BOV line- blanketed 

stellar-atmosphere model of Hickok and Morton (1968), taking the 
stellar radius to be 8 R q . It will also be assumed, following Herbig, 
that the -15 km sec ^ cloud is an H I region which is at least 15 pc 
from Q Oph, thereby setting an upper limit to the contribution of 
Q Oph to the total radiation field; 
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Consider now the attenuation of these fields by the cloud. The 

color excess of Q Oph has been determined by Bless and Savage (1970) 

to be E = +0. 28, which according to Whitford (1958) represents a 
B- V 

o 

total extinction of 0. 87 mag at A5480 A. It is reasonable to assume 
that this extinction is produced mainly by the -15 km sec ^ cloud, and • 
that the extinction at the center of the cloud is 0. 44 mag. The Q Oph 

o 

reddening curve has a hump at \2200 A and increases with inverse 
wavelength somewhat faster than the customary \ ^ law. If the grains 
within the cloud have appreciable albedo, however, scattered starlight 
will reduce the attenuation in the cloud markedly from that expected if 
the extinction were purely absorptive. Habing (1968) gives a simple 
approximate expression for the attenuation factor G as a function of 
albedo a and the total extinction A(\) in a homogeneous cloud, namely 

-1 

(3) 

New observations of the general galactic field by Witt and Lillie (1971) 
indicate that the albedo has a minimum of 10 percent at the position 

o 

of the extinction hump (A.2200 A) and rises sharply to 90 percent at 
X1500 A. In the short wavelength region, then, where the albedo is 
near unity there is not much absorption of photons even though the 
total extinction may exceed 3 magnitudes. 

Since the actual position of the -15 km sec * cloud is poorly known, 
we are led to adopt two cases which hopefully bracket the range of possible 


CT = 


a + (~l-a)10 


0. 4A(M 
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radiation densities: (1) a "high" field which is the sum of the galactic 
and the stellar radiation at 15 pc from Q Oph, attenuated by the 
factor G computed from the reddening curve of Bless, and Savage 
and the albedo data of Witt and Lillie; and (2) a "low" field which is the 
attenuated galactic radiation alone. These "high" and "low" fields 
as a function of wavelength are shown in Figure 7. 

iii ) Atomic Photoionization Cross-Sections 
and Recombination Rates 

For Ca I we follow Herbig in piecing together the photoionization 
cross-sections at different wavelengths of Hudson and Carter 
(1967a), Newsom (1966), and Ditchburn and Hudson (I960): the 
composite is shown in Figure 8. For Na I and K I, we use the 
cross-sections of Hudson and Carter (1967b and 1965). For the 
A1 I cross-section we adopt the expermental data of Kozlov, Nikonova, 

o 

and Startsev (1966) in the region AA1550 — 2070 A, and the calculation 
of Burgess, Field, and Michie (I960) at shorter wavelengths. For 
Fe I, we use the recent theoretical cross-section of Kelly and Ron 
(1971), based on many-body perturbation theory, which is consider- 

• o 

ably larger than the Hartree-Fock value in the region X912 — 1970 A. 

The recombination rates of Ca I, Na I, and K I are taken from 
Seaton (1951) and that of A1 I from Burgess, Field, and Michie (I960). 
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The recombination rate of Fe I is highly uncertain; a value 

0. 6 x 10 sec \ comparable to that of the other atoms, will be 

adopted. 

Table 4 lists the atomic photoionization rates computed from 
Equation (2) for the two adopted radiation fields and the adopted 
recombination rates. 

iv) The Electron Density 

In order to calculate the electron density with Equation (1) 
we make the additional assumption that low density matter at 
-15 km sec ^ exterior to the cloud does; not contribute to the observed 
equivalent widths of the calcium lines, an effect which Herbig estimates 
will influence the column densities by less than a factor 2. Equation (1) 

then yields with the Ca I /Ca II line ratio and the appropriate r/a an 

- 3 - 3 

electron density equal to 0. 55 cm for the high field and 0. 16 cm 

for the low. If the albedo data were grossly in error, the electron 

- 3 

density could be as high as 0.75 cm —corresponding to a completely 

- 3 

unattenuated high field — or as low as 0.07 cm —an attenuated low 
field where there is no scattered galactic light . We consider neither 
of these cases likely. In his calculation of the electron density, Herbig 
(1968) assumed essentially the unattenuated high field at a star- cloud 
distance of 20 pc. ' Using the assumed normal abundance of sodium 
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relative to hydrogen to estimate N(Na II), he obtained the electron 

-3 

density n g = 0. 3 5 cm 

The total column density of sodium may now be determined 

without Herbig's assumption that the Na/H ratio has the solar system 

value. Taking the column density he determined from the observed 

13-2 

D and D 7 lines,. N(Na I) = 5 x 10 cm , Equation (1) yields for the 

14 - 2 

total sodium column density N(Na) = 6,1 x 10 cm . For the hydro- 

20 - 2 

gen column density Stecher (1968) finds N(H I) = 4. 2 x 10 cm from 
rocket measurement of Lyman a. The Na/H abundance ratio is then 
found to be 1. 5 x 10 which is indeed consistent with the solar 
system value 2.4 x 10 ^ (Cameron 1968). 

It should be remarked that N(Na) is insensitive to the adopted 

radiation field. On substituting n = F N(Ca I)/a_ N(Ca II) into 

G C C Si 

the sodium ionization equation it becomes apparent that the only 

dependence of N(Na) on the radiation field is through the ratio 

r /r„ • As may be seen from Table 5, this ratio and also the 
IMa Ca 

ratios K/Ca, Fe/Ca, and Al/Ca are insensitive to the exact 
radiation field. 

v) Atomic Abundances 

Table 5 lists the atomic abundances computed from Equation 
( 1) and the data in Tables 2 and 4. The values quoted are an average 

25 
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of those found for the high and low fields which differ typically by 
less than 20 percent. Sodium and potassium relative to hydrogen 
have their normal solar system abundances, within the observational 
uncertainties. Gaseous calcium here, as elsewhere (cf. Howard, 
Wentzel, and McGee 1963), seems to be deficient by a large factor. 
The factor 2500 here is, however, somewhat greater than generally 
observed. ^ 

Iron and aluminum also appear to be deficient by large 

factors in the gas phase. From the Fe I column density, 

N(Fe I) = 3. 5 X 10 11 cm and the data in Table 4, the total 

13 - 2 

iron column density is found to be 4. 3 x 10 cm . The Fe/H 

-7 

ratio is then 1., 0 x 10 , almost three orders of magnitude below the 

solar system value. This discrepancy is probably not attributable to 
uncertainty in the calculation of the photoionization cross-section 
(Kelly 1971, private communication). One cannot evaluate precisely 
the uncertainty in the recombination rate, but it is presumably 
less than an order of magnitude. 

An important result of the observations is the apparently large 
underabundance of aluminum in the -15 km sec * cloud. The 0. 26 mA 

o 10 - 2 

upper limit for A1 I \3944 A yields N(A1 I) < 6 x 10 cm . We 
obtain an upper limit on the total column density of aluminum 

! _ 2S 
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N(A1) < 1. 1 x 10 cm , which then gives an Al/H ratio of less than 
- 8 

2.6 x 10 , a factor of at least 100 less than the fractional abundance 

in the solar system (Cameron 1968). 

It is unlikely that new atomic data will remove the large 
deficiencies of iron and aluminum, and our results consequently 
support the idea that most of these metals is trapped in the inter- 
stellar grains. Herbig (1968) has noted also a titanium underabundance 
in the -15 km sec ^ cloud by a factor greater than 100. The question 
of deficient calcium has not been satisfactorily answered. For the 
radiation fields adopted, the ionization of Ca II should be negligible 
and other depletion mechanisms must be considered. Goldsmith, Field, 
and Habing (1969) have attempted to explain the low abundance of calcium 
relative to sodium by differences in their rates of accretion on grains, 
but found that the differences are not sufficient to produce the required 
Ca/Na ratio without seriously depleting both these elements relative 
to hydrogen. 

vi) The Hydrogen Density and the Degree of Ionization 

It is possible to estimate the neutral hydrogen density from the 
electron density if we assume normal cosmic abundances of the 
predominantly once-ionized elements carbon, magnesium, silicon, 
and iron. Change neutrality gives 
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n = n(H II) + n(C n) + n(Mg II) + n(Si II) + N(Fe II) 
e 


(4) 


the right-hand side of which may be written n(H II) + en(H I), where 
-4 

e = 6 X 10 according to the abundances given by Cameron (1968). 
The ionization equilibrium for hydrogen is 


n(H II) n 

± _ X 

n(H II a 


(5) 


where Q is the total cosmic ray and soft X-ray ionization rate and 
a is the total recombination rate to all but the ground state. (Since 
the cloud is assumed to be an H I region, there is no photoionization 
of hydrogen .) According to Seaton (1959), a = 0.69 x 10 ^sec \ 
Therefore we arrive at an expression for n^ in terms of n(H I) and 


n 

e 


en(H I) 
2 


1 ' 
2 


2 2 1 1/2 
e n(H I) + 4£n(H I) /a 


( 6 ) 


The energetic particle ionization rate Q is highly uncertain. 

Spitzer and Tomasko (1968) give the upper and lower limits 

1. 17 x 10 ^sec * and 6.8 x 10 ^sec ^ for ionization due to 

cosmic rays, and a recent calculation by Hjellming, Gordon, and 

Gordon (1969), based on dispersion of pulsar signals, indicates that 
-15-1 

Q 2. 5 x 10 sec . Silk anc^^rner (1968) estimate that the 


C- 
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ionization rate due to X-rays which have passed entirely through 

-17 -1 

the hydrogen in the galactic plane is about 8 X 10 sec — 
probably small with respect to that of cosmic rays. 

A further uncertainty that must be considered is the possible 

depletion of the easily ionizable elements, notably carbon, by 

accretion on grains. A set of depletion factors has been computed 

by Goldsmith, Habing, and Field (1969): they find the factors 2, 10, 

5, and 3 for O I,. C II, Si II, and Fe I, respectively, at an accretion 
7 

time of 7 X 10 years. (The ions accrete more rapidly since collisions 
with fast moving electrons produce a negative charge on the grains.) We 
will assume that the Mg II is depleted by the same factor as Si II. 

The deficiencies of iron and aluminum noted in the preceding section 
are highly suggestive of this depletion. Smith (1972, to be published) 
has obtained a rocket UV spectrum of £ Oph which indicates that / 
carbon is depleted by a factor no more than five. Thus the depletion 
factors of Goldsmith et al . (1969) for the major trace elements lead 
to a reasonable lower limit on their abundance with respect to hydrogen. 

_4 

These lowered abundances lead to e = 0.7 x 10 . Although 

the effect of this decrease on n is quite large when £ is- small, 

- 15 

it is minimal up to rather high values of n(H I) when £ = 1. 17 x 10 
sec . Figure 9 shows n^ versus n(H I) for the two values of g 
and several values of £. We distinguish two cases' 
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1. Normal e • The effect of cosmic ray ionization is 

— 3 -3 

small for n(H I) > 300 cm ; for n(H I) = 10 cm , cosmic rays 
produce nearly a factor 10 more electrons than the easily 

_ 3 

ionizable elements. Given the limits n = 0. 16 to 0. 55 cm , 


e 

-3 

n(H I) falls in the range 260 —900 cm for Q = 0, 90 — 570 cm 
-15-1 -3 

for Q= 1. 17 x 10 sec , and 50 — 440 cm for Q = 2.5 x 10 


2. Depleted e. The effects of cosmic ray ionization are 

qualitatively the same. For Q = 0, n(H T lies between 2200 and 
- 3 - 15 - 1 

7900 cm ; for Q = 1. 17 x 10 sec , between 130 and 1500 cm 

- 15 - 1 -3 

and for £ = 2. 5 x 10 sec , between 60 and 740 cm 

One ..objection to a .large hydrogen '"density is the resulting 

-3 

thinness of the cloud. For a density of 1000 cm , the thickness 
is computed from the hydrogen column density to be about 0. 15 pc. 
This is to be contrasted with the fact that the -15 km sec ^ cloud 
extends over tens of parsecs perpendicular to the line of sight. 

As Habing (1969) has pointed out, a rough estimate of the lifetime 
of such a thin cloud based on the observed turbulent velocity (that 
determined from the curve of growth doppler constant) yields a 
value comparable to the time needed to establish ionization equilib 
rium. Thus, for such a thin cloud there is the possibility that any 
estimate of abundances based on the assumption of ionization 
equilibrium may be in error. 
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A more satisfying argumert for a low value of n(H I) than the 
preceding may be derived from an investigation of the energy balance 
between the cosmic ray heating and the cooling rates in the cloud. 

I 

Two different models for interstellar H I regions based on this 

approach have recently been proposed (Field, Goldsmith, and Habing 

1969 [FGH]; Hjellming, Gordon, and Gordon 1969 [HGG]). FGH 

compute the neutral particle density as a function of temperature 

- 15 - 1 

and trace element depletion for ( = 0, 6 x 10 sec . At 100° K 

-3 

they find n(H I) to be 0. 7, 9, and 32 cm , for carbon depletion by 
factors of 1 (normal abundance), 10, and 100, respectively. HGG intro- 
duce ,the hydrogen ionization equation in such a way that the ratios 

n /n(H I) and £/n(H I) can be calculated separately as functions of 
e 

temperature alone. For 100° K their model predicts n /n(H I) = 

0.01 and c/n(H I) = 0. 9 X 10 ^sec ’'’cm ^ , leading in the - 15 km sec 

-3 

cloud to n(H I) between 16 and 55 cm and to q in the rather high 

-5 -1 

range 14 to 50 x 10 sec . The HGG model seems to some degree 
inconsistent with the fairly high electron density found in this cloud; 
this inconsistency grows worse as the coolants, particularly carbon, 
are depleted. Since it is not yet possible to decide between the two 
models by observations (Hobbs and Zuckermann 1972), we will, for the 
purpose of finding an upper limit to the neutral hydrogen density, 

31 



25 


adopt the HGG results and assume that they scale with carbon depletion 

in the same manner as those of FGH. This gives an upper limit N(H I) 

-3 

= 700 cm . If the FGH model should prove the more correct, then 

_3 

the limit will be considerably less — below 100 cm , even at a 
kinetic temperature of 30° K. 

-3 

Taking the lower limit on n(H I) to be about 50 cm from the 

discussion of the hydrogen ionization and the upper limit discussed above, 

» _ 3 

we finally conclude that n(fi I) probably falls in the range 50 — • 700 cm . 

t 

vii ) Metastable Helium and the Cosmic Ray Ionization Rate 


It has been pointed out recently that there may be sufficient 

3 

interstellar helium in the 2 S metastable state to cause detectable 
absorption at A3889 A (Scherb 1968, Rees, Sciama, and Stobbs 1968). 

This state would be populated by recombination to triplet levels following 

ionization of helium by cosmic rays. According to Rees et al . (1968), 

3 8 — 3 

the life time of the 2 S state is 2. 5 X 10 /[l+2. 8 X 10 y n(H I)] sec 

where the yet unmeasured helium excitation cross-section is given by 

-16 -2 . 8 
y (y £ 5) in units of 10 cm . One thus expects the fraction 1. 8 X 10 

-3 

£ /[1+2. 8 X 10 y n(H I)] of the helium atoms to exist in the metastable 

o 

state. Hence an observation of the A 3889 A line which arises from this 
state will set limits on , which is approximately twice the cosmic ray 


ionization rate of hydrogen. Our upper limit on A 3889 A, 0. 25 mA 
(Table 1), indicates a column density N(He I*) £ 3 X 10 1( ^ cm 
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If helium is normally abundant, then N(He I) - 0. 1 N(H I) - 
4 X 10 19 cm' 2 , or N(He I*)/N(He I) <. 7 X 10" 10 and 

IQ _ 3 _ 1 

r £ 4 x 10" [ 1 + 2.8 x 10" y n(H I)] sec . (7) 

b He 

For y = 5, this gives an upper limit on which varies between 

-18 -17 -1 

7 X 10" and X 10 sec as n(H I) varies between 50 and 
-3 *‘ 

700 cm — the limits obtained above. The larger figure agrees 
with the lower limit on cosmic ray flux computed by Spitzer and 
Tomasko (1968), and therefore contradicts the result of Hjellming 
et al. (1969) based on pulsar observations, but as the lifetime of 
the metastable state is not well known, one must not take this 
discrepancy too seriously. 

d) Interstellar C^H + 

As mentioned earlier, a new interstellar line has been detected 

just to the blue of CH + A.4232. It cannot easily be ascribed to a high 

velocity component of the common isotopic species and is probably 
© 13 + 

the A.4232 A line of C H (Bortolot and Thaddeus 1969). This line 

was sought first by Wilson (1948) and more recently by Augason and 

Herbig (1967) to determine the important isotopic abundance ratio 
12 13 

a = [C ]/[C ]. The terrestrial value of this ratio is a = 89, while 

the equilibrium value in tie CNO bi-cycle is a«4 (Fowler, Caughlan, 
and Zimmerman 1967). It has been known for some time that a 
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comparably low ratio is observed in a number of low temperature 
carbon stars (McKellar 1965). From a synthesis of three spectra 
(also of Q Oph obtained with the 120 -inch coude spectrograph) Augason 

o 

and Herbig were able to set an upper limit of 1. 5 mA for the equivalent 

width of this line. This limit yielded a > 30, suggesting that 
13 

the C -rich carbon stars contributed little to the interstellar 
medium in the direction of Q Oph. ». .. . 


The isotope shift is a sum of two components. The first is 
the difference in vibrational level energies due to the change in 
molecular reduced mass, and the second is the electronic shift, 
which is generally much smaller. From Herzberg (1950, Equations 
III, 195 and III, 203), the isotope shift of a diatomic molecule may 
be written 


A\ = 


j(l-p) |V (v' + |) - t/(v" + §)] - (l - p 2 ) 


/ / , / , 1.2 // u . ii i.2 / 2 \ 

(U x (v +i?) - uj x (v +g) + 1 1 -p I 

e e e e ' 


B* J'M' + D-B" J" (/ +1) 
v v 


( 8 ) 


where p = 0. 994022 is the ratio of the normal to the isotopic re- 
duced mass, and cubic and higher terms and the electronic shift 
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are neglected. Following the usual convention, ' refers to the 

upper and " to the lower state. Herzberg (1950) gives 

6 ' 

uu" - 2((jj" x) = 2739. 54 cm" 1 
e e e 

s' 

..uu ' = 1850. 02 

e 

*' uu ' x = 103. 93 

e e 


and 


B = 11. 43 

o 


The quantity a) x is unknown for CH , but we will adopt the value 
e e 

-1 U. 

49 cm , which is uu x of B H and must be close to the actual 


e e 

value. Thus, from Equation (8), we have AA. = - 0. 255 A. The 
more recent molecular constants of Douglas and Morton (I960) 
yield the same result. 


The electronic shift between C^h"*" and is small. 

According to Douglas (1972, private communication), it should 
be less than 15 percent of the isotope shift computed above. He 
also estimated that the errors in calculating the usual isotope 
shift (that is, Equation (8)) are about as large as the electronic 
shift. The sum of these smaller effects, however, amounts to 
a total error of less jthan one line width and may be neglected. 
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Figure 3 shows the results of the present synthesis in the 

o 12 + 

vicinity of the \4232 A line. The strong C H line falls at =* 15. 0 
km sec . Taking the isotope, shift computed above, the new line 
falls at -15. 2 km sec ^ with an equivalent width equal to 0.4 ± 

o 

0. 13 mA. . 

The signal strength is over three times the standard deviation 

of the noise in the filtered spectrum, so there is little doubt that 

the feature is real. It is hard to see how it could result from an 

adjacency effect, and in any case no spurious line of this nature 

has been observed in the vicinity of the other strong lines. We 

must, however, consider seriously the possibility that the new 

. 12 + 

line is simply a high velocity component of the C H line. It has 
been noted that the H and K lines are double at high dispersion, 
with a high velocity component at -29 km sec \ This velocity 

o’ 12 V 

corresponds to a -0. 19 A displacement from the C H line, more 

than a linewidth short of the new feature. Thus it is implausible 

13 + 

that the line which we attribute to C H is produced wholly in the 
-29 km sec 1 cloud, although the possibility that it is blended with 
a line at this velocity cannot be entirely discounted. Hobbs (1964, 
private communication), with his Pepsios spectrometer has observed 
a weak component of sodium at -31. 6 ± 0. 2 km sec . A molec- 
ular cloud of this velocity would produce a line at -0. 25 A with 
respect to C^H + , virtually at the wavelength of the new feature. 
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We reject this idea, however, on the grounds that it would 
require an interstellar CH + /Na density ratio far in excess of any 
so far observed. Both the new line and Hobb' s -31.6 km sec ^ 
sodium line (W « 1 ml) are undoubtedly unsaturated, and the ratio 

of equivalent widths W(CH + )/W(Na I) «s 1/2 is therefore the ratio 
of optical depths. On the other hand, in the -15 km sec * cloud this 

ratio is 1/65 (Herbig 1968), typical for those clouds where 14232 A 
is conspicuous, although there are admittedly few clouds where 
the curve of growth is well enough known to permit this ratio to be 
calculated precisely. However, by taking a reasonable range of 
values for the doppler constant b, say 1 to 4 km sec \ it is possible 
to extract approximate values of t(CH + )/-r(Na I) from available data. 
Thus, for 34 of the 6l stars for which W(CH + ) has been measured 
or estimated by Adams (1947) or by Rogerson, Spitzer, and Bahng 

(1959), W(Na I) is also found in the various compilations of inter- 
stellar Na I and Ca II lines (Merrill et al. 1937; Spitzer 1948; Beals 

and Oke 1953; Munch 1957). Among these, there are a dozen cases 
where the structure of the D, or H and K, lines suggests that the 
absorption takes place in essentially a single cloud, in which the 
CH + line is presumed formed, on the basis of common radial velo- 
cities. In all these cases, W(CH + )/W(Na I) < 1/5, while if the 
high degree of saturation in the line is taken into consideration, 
the optical depth ratio becomes markedly smaller, r(CH + )/T(Na I) 
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< 1/17, assuming b in the range 1 to 4 km sec . While further 
inves tigation of this ratio is indicated, the data just cited are suffi- 

o 

ciently extensive to suggest strongly that A4232.08 A is unlikely to 

12 + 

be a high velocity component of the C H line, and we instead 

assign this new feature, on the basis of its precise location at the 

13 + 

expected isotope shift, to the species C H . 

13 + 

A potentially conclusive test for the existence of C H in 
the -15 km sec ^ cloud would be detection of the equivalent line at 
A3958 A, which is the CH + R(0) line of the(l, 0) band in the same 

O 

electronic system. The isotope shift of this line is +0.44 A and 
there can be no confusion with known high velocity clouds. The 
Franck-Condon factor for this band is favorable, being roughly 
half that of the A.4232 A band according to Herbig (1968). No such 
line is found, but the upper limit at the wavelength indicated is 

0. 25 mA., consistent with the idea that such a weak line as expected 
is just below the threshold of detectability. The same situation 
holds for C^N. 

The optical depths of the C 12 H + and C 13 H + lines are t = 

1 2 

0. 968 and = 0* 012 ±0. 003 from the curve of growth. If we 
rather pessimistically assume that a third of the C* 3 H + line may 
be due to a blend with the stronger line, then r = 0. 012 (+0. 003, 
-0. 007) and the abundance ratio in the line of sight is 
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a = t 12 / t 13 = 81 (+113, -16) . (8) 

13 

Recent observations of interstellar CO and C O and Q Oph have 
yielded as comparable value of a (Smith and Stecher 1971). Our 

finding that a is essentially terrestrial in the - 15 km sec * cloud 

13 + 

clearly suggests that C -rich stars have not contributed to the CH 

molecules, or presumably to the rest of the interstellar matter, . 

in the direction of Q Oph, and moreover it implies that the 

composition of the interstellar medium has remained rather 

9 

stable for the age of the solar system, or 5 X 10 years. 
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III. THE + 13 KM/ SEC" 1 CLOUD IN Q PERSEI 

The Perseus II association of which Q Per is a member has 
been studied in considerable detail (see Lynds ( 1969 ) for references), 
but the interstellar line spectrum of £ Per, though extensive at 
+ 13 km sec \ has received relatively little attention until re- 
cently. The high resolution Fabry-Perot work of Hobbs (1969 and 
1972, private communication) has resulted in a fairly good idea 
of the curve of growth and White (1972) has possibly detected 
Ca I \4227 A with the coude scanner of the 100-inch telescope at 
Mt. Wilson. Since the position of the + 13 km sec ^ cloud relative 
to Q Per and the other stars in the association is uncertain, 
it is difficult to establish the ionization equilibrium. Even so we 
may to some extent proceed as we did for Q Oph and make an 
attempt to determine some of the physical conditions and abundances 
in the cloud. 

a) Observational Results 

Seven high dispersion spectrograms, similar to the Q Oph 
plates described earlier, were obtained with the 120-inch telescope 
and its coude' spectrograph. .The long wavelength half of four of these 
was not in proper focus due to a maladjustment in the spectrograph, 
and were not analyzed. The good spectrograms were reduced 

40 


34 


to intensity versus wavelength and synthesized by the method 
described above and in Appendix A. 

The observed atomic and molecular line profiles are shown 

in Figures 10 and 11, and the. vicinity of the CN lines in Figure 12. 

Three of the atomic lines are seen in C Per for the first time. These 

are Fe I X3720 and X3860 A and Ga I X4227 A, And all fall at + 13 

km sec , the velocity of the strong interstellar lines. Although 

Ca I X4227 A is barely detected, it is noteworthy that the line is 

.also observed. by White (1972), with the same equivalent width. 

The strong H and K lines of Ca II appear triple at this dispersion 

with the main component at + 13 km sec ^ and weaker lines at 

-1 = , 

-6 and + 21 km sec . The A 3303 A D lines of Na I show similar 
structure on four 120 -inch coude spectrograms kindly provided 
by Dr. L. Kuhi. Table 6 lists the equivalent widths and optical 
depths of all lines detected in the spectral interval XX 3650 -4350 A. 
This table includes as well upper limits to lines in the CH, CH, and 
CH + bands, of interest in the measurement of the microwave back- 

o 

ground radiation, and to the A1 I X3944 A line from which an upper 
limit to the aluminum abundance will be derived. 

As in the case of Q Oph there appears to be no systematic 
difference between atomic and molecular radial velocities. The 

' ' .. , 41 : . ... , 
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slight inaccuracy in the CH + A3958 A. laboratory wavelength 
suspected in Q Oph seems to be confirmed in this cloud. 

The curve of growth in the + 13 km sec ^ cloud determined 
from the synthesis alone is quite uncertain. The doublet ratio 

method yields b(Ca II) = 2. 1 km sec" 1 for the H and K 

alines, while b(Na I) calculated from the D 7 lines is 0. 7 km sec \ 

This discrepancy is only partially explainable by errors in the 
equivalent widths. The sodium doppler constant may also be 

obtained from the ratio of the to the D 7 ^ line (Stromgren 1948); 

. -1 

Hobbs (1969) gives W (D^) ^ rn ^ > so that b(Na I) =1.6 km sec 

Unfortunately, the relatively poor agreement between photographic 

and photoelectric equivalent widths makes the second calculation 

uncertain. However, the lower value, 0. 7 km sec \ is found also 

o 

for the CH doppler constant computed for the A. 3886 and A 4300 A lines 
using the electronic f- values of Fink and Welge (1967) and the Franck- 
Condon factors given by Herbig (1968). The CH doppler constant is 
strongly affected by uncertainties in the f-values and may be as high 
as 1.5 km sec ^ . Hobbs (1972, private communication) has, however, 
resolved CH + A4232 and reports its doppler constant to be 1. 98 km sec 
corrected for instrumental profile and saturation. We will adopt this 
value for all the molecules. To sum up, the atomic and molecular 
column densities (Tables 7 and 8) have been computed from the curve 
of growth with the doppler constants 2. 1 km sec ^ for Ca II, 0. 7 
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km sec for Na I, and 1. 98 for the molecules. For the faint lines, 
no curve of growth correction is necessary. 

b) The Electron Density and the Atomic Abundances 

Owing to a lack of very definite knowledge about the location of 
the + 13 km sec ^ cloud with respect to f Per, it is difficult to specify 
the radiation field at the cloud center. There are some grounds, 
however, for placing the cloud at a considerable distance from Q Per. 
The following arguments are based primarily on the observational data 
summarized by Lynds (1969). 

Galaxy counts in the region of the Perseus II association (Shane 
and Wirtanen 1967) indicate three regions of different obscuration. 

Of the six Perseus II stars in Adam's (1949) survey which show molec- 
ular absorption, three — Q Per, X Per, and o Per — fall in the region 
of highest extinction. These exhibit a significantly different CH + /CH 
line ratio from the three stars which fall outside — § Per, HD 22951, 
and HD 24131. It is reasonable to associate the molecules observed 
in Q Per and the other two stars with this high extinction region, and on 
the basis of their radial velocity, it is reasonable to assume that the 
atomic lines are formed there as well. According £o Lynds (1969), 
the best estimate of the distances of the heavily obscuring material 
is 200— 300 pc from the sun while the association is believed to lie 
some 300 to 400 pc away. Thus, the cloud and the stellar association 
which contains q Per may be 100 pc or more apart. 
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Of the three stars lying behind the obscuring layer, the nearest 
is 6 Per. The distance moduli determined by Blaauw and Borgman 
(1963), 7. 59 mag for Q Per and 6. 48 mag for o Per, correspond to 
distances of 340 and 190 pc. According to Stothers (1972 private com- 
munication) Q Per is a peculiar supergiant and its absolute visual 
magnitude may be in error by as much as 0. 5 mag, but even so, the 
large difference in distance moduli between the two stars retains 
some significance. Most likely then, the nearest hot star to the cloud 
under study is o Per (MK type B1III). If o Per is 190 pc from the 
sun, then it is at least 13 pc from the column seen against Q Per. 

We may now proceed to estimate the radiation field at the cloud 
center. As before, we adopt as a lower limit on the radiation field 
the diffuse galactic field proposed by Habing (1968), attenuated by 
the cloud. As an upper limit we take the sum of this field and that of 
0 Per at 13 pc, attenuated to a lesser extent since we have really no 
idea of the extinction between this star and the cloud center. For 

o 

the total extinction we adopt 0. 54 mag at A5480 A, half that measured 
by Bless and Savage (1970), and combine their £ Per redden- 
ing curve (that of q Oph is used at wavelengths above their region 
of measurement) with the Witt and Lillie (1971) albedo data to 
obtain the attenuation factor. To compute the stellar field 
attenuation we will adopt the smaller extinction value 0. 27 mag 
since the extinction between the cloud center and 0 Per may be 
smaller than that along .the line jaf sight. The line -blanketed B2II-III 
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model of Van Citters and Morton (1970) provides suitable fluxes for 
calculation of the stellar field. The stellar radius is taken to be 
13 R .' The atomic photoionization rates for the two fields adopted 
are given in Table 9. 

Assuming the existence of Ca I 14227 A, the Ca I/Ca II line 

-3 

ratio yields an electron density in the range 0. 044 to 0. 14 cm 

Further observations are obviously needed to substantiate this result 

The neutral hydrogen densities indicated by this range of n may 

- 15 - 1 

be read from Figure 5. For Q = 1. 17 x 10 sec , n(H I) lies 

-3 

between 10 and 80 cm , if carbon is depleted by a factor no 
greater than 10. 

The total column densities of the atoms found in the 
+ 13 km sec ^ cloud are given in Table 10. Since all the elements 
are presumed to be predominately once-ionized, their relative 
abundances will probably not change greatly if the electron density 
found above is in error. It is of some interest to compare the rela- 
tive abundances observed here and in Q Oph. N(H I) has not yet 
been measured against Q Per, so the abundances shown in Table 11 
are defined relative to sodium whose fractional abundance relative, 
to hydrogen is known to show little variation in the interstellar 
medium. The deficiencies in iron, calcium, and aluminum are 
found to be similar in the two clouds. 
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IV. THE SPECTRUM OF THE MICROWAVE 
BACKGROUND RADIATION 

The history of the microwave background radiation began 
in 1948 when Gamow and his coworkers, in the course of their 
investigation into cosmological element synthesis, first recognized 
the importance of hot thermal radiation in the early expanding 
universe . (Alpher, Bethe, and Gamow 1948; Gamow 1948). They 
predicted that a highly red- shifted remnant of this radiation would 
survive into the present era with its blackbody spectrum intact, 
and made estimates of its temperature falling in the range 
5 — 28° K (Alpher and Herman 1949, 1951). Although it is 
conceivable that this radiation could have been detected with 
the radio receivers of the day, no deliberate search was attempted 
until nearly two decades later. 

At the suggestion of Dicke in 1965, Roll and Wilkinson at 
Princeton started to construct a 3. 2 cm receiver for detection 
of Gamow' s radiation. In the meantime, while calibrating the 
20-foot horn antenna of the Bell Telephone Laboratory at Holmdel, 
New Jersey, Penzias and Wilson (1965) discovered an excess 
antenna noise at X7.35 cm which corresponded to a temperature 
of about 3. 5° K. This noise, to their experimental accuracy, 
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was isotropic and unpolarized, and "showed no time dependence. 

Their discovery and its interpretation were published in the spring 
of 1965 (Penzias and Wilson 1965; Dicke, Peebles, Roll and 
Wilkinson 1965). 

The first confirmatory evidence that the radiation had the 
predicted blackbody character came with the realization that the 
classical CN observations of Adams (1941) implied ~3° K radia- 
tion at \ 2 . 6 mm. This was rapidly followed by Roll and Wilkinson's 
(1966) \ 3 . 2 cm measurement, which yielded a brightness tempera- 
ture of 3. 0 ± 0. 5 °K, and a series of further direct measurements 
that established the blackbody character (T » 2.7° K) of the 
radiation over more than two decades of wavelength, from \73. 5 cm 
to \3. 3 mm (see Thaddeus 1972 for references). 

A crucial region of the spectrum is near 1 mm wavelength, 
where ~ 3° K radiation has its peak and most of the energy lies. 

This region is inaccessible from the ground because of the opacity 
of the atmosphere, but in the past four years there have been direct 
observations with rockets and balloons from above the atmosphere. 

These observations disagree among themselves to a great extent 
and have led to considerable controversy concerning the intensity 
and shape of the background spectrum. The measurements from 
the interstellar molecular bands, which concern us here, are 
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conservative in nature, and are not dependent on the accuracy of 
large instrumental and atmospheric corrections. Although the effect 
of the curve of growth correction on the molecular upper limits is in 
fact sub stantial, it is rather accurately known and any error in it 
is more likely to still further lower the upper limits. The molecular 
bands provide one precise measurement and three fundamental con- 
straints at millimeter and submillimeter wavelengths against which 
all direct observations must be compared. 

Much of the observational material on CN, CH, and CH + in 
Q Oph described here was analyzed in a preliminary way by Bortolot 
.et.al. (1969). S.ome additional plates have since been obtained, and the 
densitometry is completely new and more refined. The analysis has 
further benefitted from a better curve of growth and a more detailed 
statistical treatment of upper limits to missing lines. 

If the CN excitation observed is indeed due to a universal back- 
radiation, it must show directional invariance. The data available 
regarding this will be discussed at the conclusion of this section. 

a) The Molecular Excitation Temperatures 

The excitation temperature of a molecular band may be de- 
termined directly from the intensities of any two absorption lines 
that do not share a common lower level, providing the populations 
can indeed be described by a single temperature. Consider two 
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optical absorption lines which arise from transitions from electronic 
ground state levels i and j (E. > E.) to excited electronic levels i 7 
and j 7 . The transition intensities will be proportional to n, and 
n., the populations of levels i and j. Alternatively, the population 
ratio may be inferred from the observed line strengths, being given by 



g.S.., T.. / 

1 JJ 11 

g.S.., T../ 
J ii JJ 



where the g's are the statistical weights of the levels, the 

S's the theoretical line strengths (the summed squared transition 

matrix elements), and the t's the optical depths of the lines. If 

we define s = g.S..//gS../ and r = T ../ / t. J , then n./n. = r/ s 
J ii i JJ 11 JJ i J 

and the excitation temperature is given by 


T = hv /k^(s/r) (10) 

ij ij 

where h and k are the Planck and Boltzmann constants^ and 

hv.. = E. - E. . Figure 13 shows partial electronic and rotational — 
ij i J 

or in the case of CH, fine structure — term diagrams of the three 
molecules observed optically in interstellar space. The wave- 
lengths at which we will attempt to determine the brightness 
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temperature of the background radiation are X2. 64 and XI. 32 mm 
for CN, and XO. 56 and XO. 36 mm for CH and CH + , respectively. 

Direct measurements of the background radiation from high 
altitude sites may ultimately be made in the XXI. 6-2. 6 mm radio 
window, but these will undoubtedly be extremely difficult. There 
does not appear to be any prospect of direct measurements from 
the ground at X £ 1 mm. 

The requisite line strength ratios s are obtained for CN and 
CH + from the Honl- London factors (see Herzberg 1950, p. 208): 


J. + J.' + 1 


J. + J.' + 1 
J J 


( 11 ) 


where the J's are the total angular momenta of the levels indicated 
by subscripts. The CH + ground state is ^2; that of CN is but 

the doublet splitting is unresolved optically, and Equation (11) 
is valid. In the case of CH, which is intermediate between 
Hund's cases (a) and (b), the calculation is somewhat more 
involved. It requires a transformation of the matrix elements 
of the electric dipole moment operator, calculated in the 
pure Hund's case (b) representation, to the intermediate 
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representation where the molecular Hamiltonian matrix is diagonai. 

In Hund's case (b) one finds s = 5/3, while the correct transformed 
value is s = 1. 524. 

The likelihood that any pumping mechanism for CN or CH + 

exists which can lower T. . below T„(v. .) would appear slight for 

ij - B ij 

several reasons; the rotational level schemes are simple, the 

radiative lifetimes of the levels are short, and the density of the 

3 -3 

molecular clouds observed optically is probably not great (s 10 cm ). 

The fine structure selection rules of CH might permit pumping to occur 

under certain conditions. However, in the clouds where CH is observed 

the departure of the population ratio from that expected when the pump- 

5 

ing is ignored is estimated to be no more than one part in 10 
(Clauser 1970). Thus local processes such as collisional impact 
with charged and neutral particles or resonance scattering of 
starlight can serve only to increase the molecular excitation. 

The excitation temperatures, therefore, are almost certainly 
reliable as upper limits to the background brightness tempera- 
ture. 

The observational data on the molecules in the Q Oph -15 km sec 
cloud have been given in § II b and are shown again for convenience in 
Table 12. The errors in the CN equivalent widths represent one 
standard deviation. The optical depths were computed using the CN 
doppler constant, b = 1. 36 km sec \ measured by Hegyi et al. (1972). 
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The CN excitation temperatures, as calculated with Equation 
(10), are given in Table 13. If we make the minimum curve of 
growth correction (the observed R(0) line width in the synthesis 
sets the upper limit 2. 40 km sec ^ on the doppler constant), the 
resulting T ^ (CN) from the R(1)/R(0) and P(1)/R(0) line ratios, 

2. 99 °K, is a strict upper limit to the brightness temperature 
at X2. 64 mm. The proper curve of growth correction reduces 
this temperature by a small amount. A correction for the local 
processes in the -15 km sec ^ cloud , also fairly small, converts 
T (CN) to T (A. 2. 64 mm), as discussed immediately below. 

It is quite difficult to assign an uncertainty to the curve of growth 
correction, but reasonable limits on the doppler constant indicate 
that this uncertainty is no larger than those which result from 
statistical errors in the equivalent width measurements or from 
the collisional correction. For this reason no uncertainty is 
assigned. The final weighted mean 


T (X2. 64 mm) = 2. 75 ± 0. 10° K 

Jd 


( 12 ) 


may well be the most accurate determination of the background 
radiation temperature wavelength at any wavelength. 


Table 14 gives the upper limits to the background temperature 
and intensity at XXI. 32, 0. 56, and 0. 36 mm derived from the upper 
limits in Table 12 and these limits are plotted in Figures 14 and 15. 
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b) CN Excitation by Local Processes 


2 

Although the ground state of CN is E, the small p doubling 
of the levels, resolved in the \2. 6 mm radio lines (Jefferts, 
Fenzias, and Wilson 1970), is too small to be resolved optically, 
and therefore on the basis of spectroscopic stability, we may 
proceed as if the ground state were "*"E and the level structure 

that of a simple rigid rotator, as shown in Figure 13. Only the 
two lowest rotational levels are observed to be populated in inter- 
stellar space, so that the fractional level populations n^, n^ in 
statistical equilibrium are determined from the rate equation 




A i) + n o c oi = ° 


(13) 


where is the spontaneous radiative decay rate of the J = .1 
level 

4 2 3 
64tt V v 

A = = 1. 19 X 10" sec" , (14) 

9hc 

18 

(a = 1. 45 x 10 esu is the CN dipole moment, and = 

1. 13 x 10 11 Hz. C ^ is the radiative or collisional deexcitation 
rate from the J = 1 directly to the J = 0 level, while C ^ is the 
excitation rate from the J = 0 level to all others since the radia- 
tive cascade from higher levels to J = 1 occurs quite rapidly. 
Since the transition cross-sections for AJ > 1 in neutral particle 
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collisions will be comparable to those for AJ = 1. and C ^ 

will not satisfy detailed balance, 


'01 


10 


= ir ex p(- hv io /kT k) • 


(15) 


where g = 1, g = 3, and T is the kinetic temperature. Nonetheless, 

U X K 

while the neutral particle excitation rate may only be estimated, 
those due to electron and ion collisions are quite accurately known 
and do obey the radiative selection rules, AJ = ±1, to a good approxi- 
mation. Thus in the following treatment of collisional excitation it is 
permissible to proceed as if Equation (15) were valid. Since in the 
absence of local processes 


V n 0 = ( g l /g o) exp (- hv 10 /kT B) 


= 2. 8° K, 

j3 


( 16 ) 


Equation (13) may be solved for the excitation rate. Thus 


C 01 = ( g i /g o) A 1 exp (- hv io /kT B ) f 1 - exp (- hv 10 /kT B )] 


. , , rt -6 -1 

= 6. 1 x 10 sec 


- 1 


(17) 


This is the rate which is necessary to account for the observed CN 
excitation. 


One may also determine the correction necessary to convert 
the observed T,„ to T„ if the excitation rate from local interstellar 
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processes is known. Since exp(-hv in /kT ) will always be very 

close to unity we may set C^^/C = ^ anc * Equations (13) and (16) 

(where T replaces T ) then yield 
10 B 


3A 1 + C 01 

T 10 = hv io /k ^ — — 


(18) 


If R is the excess of over 6. 1 x 10 ^ sec \ this gives 


A T IQ = -2.0 x 10 R Q1 °K , 


(19) 


where AT is the temperature correction computed for ^3 °K 

- 6 - 1 

and R ^ is on the order of 6 x 10 sec or smaller. 


We will now consider the rates expected in the CN clouds 
which can contribute to the CN excitation, first to see whether 
any process besides photoexcitation by background radiation can 
account for the observed excitation, and secondly to determine 
what correction to T in (CN) must be made to obtain T (\2. 64 mm). 
The processes most likely to be of consequence are (1) collision 
of CN with neutrals, primarily H I, but possibly H and He; (2) 

br 

collisions with electrons; (3) collision with protons, if the CN happens 
to lie on the edge of an H II region; and (4) optical photoexcitation by 
starlight. Vibrational excitation of CN either by photons or collisions 
is too slow to be significant and need not be considered. 
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i) Neutral Particle Collisions 

While the excitation cross-sections for neutral particles on 

CN are not well known at present, they cannot greatly exceed the 

-15 2 

geometrical cross-section, which we take to be ~2 x 10 cm . 

In the - 15 km sec ^ cloud, we have concluded that n(H I) is less than 

-3 

300 cm . Assuming an upper limit on projectile velocity v = 

5 - 1 

1.2 X 10 cm sec (for 100° K H- atoms), therefore, one obtains 

C Q1 = n(H I) vct ^ 1. 6 x 10~ 7 sec" 1 , (20) 

a rate smaller than that required by a factor of about 40. If the 
effects of He I, taken at normal cosmic abundance, and H , which, as 

Ct 

Carruthers (1970) has suggested, is perhaps as abundant as H I, are 
included, C ^ increases by a factor less than 2. Hence almost cer- 
tainly, neutral particle impact is not the cause of the CN excitation. 

If there were a sufficient H I density to account for T , then 
the CO excitation temperature (Smith and Stecher 1971) and prob- 
ably the optical extinction would be much greater than observed. 

ii) Electron Impact 

The long range of the Coulomb interaction makes electron 

impact a highly effective process for rotationally exciting polar 

5 

molecules; the rate coefficient is a factor of roughly 10 greater 
than that for neutral particles at H I region kinetic velocities. 
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Since this process comes closest to accounting for the 2. 8 °K ex- 
citation, it is fortunate that both electron density in this cloud 
and the rate coefficient are rather accurately known. As we have men- 
tioned above, n in the Q Oph - 15 km sec * cloud is between 0. 07 and 

-3 

0. 23 cm . The rate coefficient has been computed by Allison and 

Dalgarno (1971) from a close- coupling calculation to be r^ = 2. 1 X 
6—3—1 

10 cm sec for 100° K electrons and at this temperature is quite 
temperature insensitive. Thus 

3. 4 x 10' 7 <: C ni = n r £ 12. 2 x 10" ? sec" 1 , (21) 

01 e 0 1 

at its maximum only a factor 6 short of the observed rate. While it 
is thus unlikely that electron impact is able to produce the observed 
excitation, this process does raise T.^ by 0. 07 to 0. 24° K. 

iii) Proton Impact 

Although rotational excitation by protons in low temperature 

H I regions is completely negligible, it is considerable in H II regions: 

Thaddeus (1972) has determined from a close- coupling calculation 

4 -3 

that a quite moderate density of 10 °K protons, 1.6 cm , is sufficient 
to account for the excitation observed. It becomes necessary, there- 
fore, to show that the molecular lines are formed outside H II regions 
if we are to discard proton impact as the excitation source. As noted 
in section II, Herbig (1968) has argued convincingly that the -15 km 
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sec cloud must be an H I region. It is probably safe to con- 
clude that proton collisions, while potentially an effective mech- 
anism for CN excitation, actually do not contribute much to it. 


iv) Optical Photoexcitation 


Photoexcitation by the general galactic field is of interest 
as a possible explanation of the directional invariance of T^(CN) 
which we have noted. The rate is greatest for the violet system, 
which _ 


C 

v 


Tte 


mhc 


V f 


for 


( 22 ) 


(Thaddeus and Clauser 1966) where u is radiation density, f = 0. 027 

- . ' A. 

is the oscillator strength of the electronic transition (Bennett and 

Dalby 1962) and \ its wavelength. Zimmermann (1968) gives the 

-18 -3« -1 -1 „ o 

galactic field u « 57 X 10 ergs cm A sec at A. 3900 A which 

is much higher than that assumed at the -15 km sec cloud center, 

hence C s 8 X 10 ^ sec \ It is obvious that this process has 

v 

virtually no effect on the CN excitation. 

Of course, light from a nearby star, in particular the OB 
star against which the molecules are typically observed, could 
increase this rate drastically, but there would then be no reason 
to expect the CN excitation to be invariant, as it is observed to be, 
and furthermore the intense UV of such stars would tend to dissociate 
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the CN. It would then be difficult to explain the existence of the 
molecule in interstellar space. 

We therefore conclude that no local process in the -15 km sec 

cloud is able to produce the observed 2. 8 °K CN excitation and that 

the correction to convert T (CN) to T (A.2. 64 mm) is in fact small 

10 15 

(0. 15 ± 0. 09 °K). The correction necessary in the Q Per +13 km sec 
cloud is smaller still (0. 04 ± 0.02 °K). 

£) Discussion of the Short Wavelength Observations 

Shivanandan (NRL,) and Houck and Harwit (Cornell; 1968) made 
the first attempt to detect the short wavelength background radiation 
directly froni above the atmosphere. Their rocket- borne far-infrared 1 
telescope picked up a large flux, apparently isotropic, about 40 times 
greater (han that of a 2. 8° K blackbody in the interval A. 1. 3 to 0.4 mm. 
Two subsequent flights by the Cornell group appeared to confirm this 
result (Houck and Harwit 1969; Pipher et al . 1971). The mean intensity 
over the instrumental passband, 5 X 10 ^ erg cm ^ sec ^ster "''Hz \ 
is considerably in excess of the molecular limits given here, 
especially that derived from CH (see Figures 14 and 15). Thus 
a radiation continuum is ruled out, and if the radiation is not 
spurious, it must be concentrated in some linelike feature or 
features which happen to avoid the molecular resonances. This 
possibility has prompted other workers to fly experiments with 
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some spectral resolution for a more detailed look at the excess 
radiation. 

The first such experiment, the balloon-borne radiometer of 
Muehlner and Weiss (1970), did provide some evidence of a strong 
line between 10 and 12 cm \ but all subsequent work has contradic- 
ted this result (Mather et al. 1-971; Blair et al. 1971). Muelhner and 
Weiss (1972) have themselves repeated their 1970 experiment with 
an improved version of their radiometer and report no line. 

Although corrected for a quite sizable atmospheric contribution, 
these new observations of Muehlner and Weiss cast considerable doubt 
on the existence of any short wavelength background radiation in . 
excess of~3 °K. The new intensities and upper limits are marked M 
in Figures 14 and 15. The results are all consistent with~3 °K 
radiation and disagree considerably with the Cornell-NR.L measure- 
ments. If the new data are taken conservatively as only upper limits, 
and if the large excess flux seen only in the Cornell-NRLi flights is 
spurious, then one is forced to recognize that there is really only 
meager evidence for the existence of any short wavelength back- 
ground radiation. 

The molecular upper limits, essentially because, as we have 
argued, neglected processes can only depress them further, provide 
the most reliable data on the intensity of the background radiation. 

It is certainly worth pursuing this direction further. The 
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shortest wavelength at which the background has been detected 
with certainty.is \ 2. 64 mm from CN, and it would be highly 
desirable to extend detection to shorter wavelengths. In principle, 
the technique of plate synthesis could be extended to permit the 
measurement of the CN R(2) line and hence the intensity at \1. 32 mm. 

Our upper limit is just a factor 2. 5 greater than the expected equivalent 
width. However, the observing time would be large and a more appro- 
priate approach would probably be to make use of photoelectric techniques. 
Hegyi, Traub, and Carleton (1972) have made the first such attempt 
to detect the R(2) line in Q Oph, using a Fabry-Perot spectrometer on 
a 60-inch telescope. If their result is interpreted in the same manner 

o 

as the present work, their upper limit is still in excess of the 0. 17 mA 
found here, but it is estimated that with a large telescope a week or two 
of observing time might suffice for the detection of this line. 

d) Invariance of the CN Excitation 

Since pumping mechanisms are excluded in the two clouds we have 
studied in detail, and we have no reason to expect pumping in the other 
clouds observed, the CN excitation temperature found in any of these 
clouds must be equal to 2. 8 °K, if its source is correctly identified with 
photoexcitation by the microwave background radiation. Such observations 
provide then a critical test of the hypothesis that the background radiation 

causes the excitation, and are the primary motive for seeking interstellar 
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CN in new stars. The existence of the short wavelength back- 
ground probably collapses if T^q(CN) can be found in a single 
case to be zero, or even appreciably less than 2. 8° K. 

So far the number of stars against which CN has been 
observed is small: Adams (1949) in his survey of 300 OB stars 
found seven, and subsequent workers have reported relatively few 
more. CN temperatures observed against several stars have 
been noted by Claus er and Thaddeus (1972) and Peimbert (1968), 
but only Q Oph and Q Per are sufficiently bright and behind 
enough CN to provide a fairly precise temperature determination. 

In this section we will refine some of the past measurements and ^ 
present the results of a search for CN against a number of late 
B and early A stars. 

£) Observational Results 

We will first discuss criteria for our choice of objects to 
observe. Experience has shown that the best place to find CN 
is against bright OB stars near the galactic plane where most 
of the interstellar material is concentrated. They have convenient 
spectroscopic properties, being intrinsically bright and having broad, 
uncomplicated stellar features as well, so most of the CN tempera- 
tures have been measured against them. These stars, however, 
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suffer from two disadvantages: those whose spectra contain CN 

are few, and more importantly, tend to be associated with H II 

regions. If the CN were located, at least in part, at the edge of 

th ionization region, then the excitation, as noted above, might be 

4 

affected by collision with thermal (10 °K) protons. If the CN 

lines are not formed unambiguously in H I regions, then the tem- 
peratures may be in doubt. Hence stars later than OB, which are 
free from extensive H II regions, would be preferable as objects 
for observation. Peimbert (1968) has discovered CN in two or 
three out of nine later-type stars he investigated. 

W.e have chosen, therefore, two groups of objects for 
study, (1) certain OB stars known to contain CN, and (2) late B 
and early A giants with extinction of ~0. 3 mag or greater. High 
dispersion spectrograms of the stars listed in Table 15 were 
obtained with the 120 -inch coude spectrograph and reduced as 
before. Of the seven later spectra, one, that of 13 Ceph, had 
the strongest CN R(0) line yet seen, excepting two distant objects 
reported by Munch (1964), and another, that of HR 7573, was 
found to contain CN weakly. 
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Except in the cases of Q Oph and Q Per, the CN doppler 
constants are not known and it is necessary to assume a value 
of b to compute optical depths. In lieu of more precise data 
we will adopt b = 1. 67 km sec \ the average value for the 
Q Oph and Q Per clouds. We will not attempt to assign an error 
due to uncertainty in the curve of growth, but since the CN lines 
are fairly weak, the curve of growth errors cannot be large 
compared to uncertainties in local excitation processes and in 
the equivalent widths themselves. The equivalent widths and 
corresponding optical depths are given in Table 16. 

The rotational excitation temperatures are presented in 
Table 17 for all stars in Table 15 in which CN is found, except 
HR 7573, where the R(l) line is not seen with any degree of 
confidence (T ^ < 3. 3° K). This table includes some additional 
objects investigated by Clauser and Thaddeus (1972) and Peimbert 
(1968). BD+66° 1674 and BD+66 0 1675 have been assigned the large 
doppler constant 5 km sec ^ by Clauser and Thaddeus on the 
basis of their great distance and low galactic latitude. Note 
should be made that the entries for Q Per and X Per represent 
measurements in essentially the same cloud, as pointed out 


in the previous section. 
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ii) Discussion 

The CN temperatures in these widely separated clouds though 
less certain than in Q Oph or Q Per are all consistent with 2. 8° K. 
This apparent invariance supports the idea of a universal back- 
ground radiation, but in addition provides what is perhaps the 
strongest evidence that the molecular lines originate in H I regions 
and that the excitation mechanism is not "local". 

While the effect of fast protons on CN excitation is exceed- 
ingly large, the lack of large variations in temperature observed 
against the early stars, when taken in the context of the great 
inhomogeneity of the interstellar medium would appear good evidence 
that the molecules do not exist in ionized regions. This idea is 
supported by the radial velocity coincidence generally observed 
between the CN and the 21 -cm lines. It is noteworthy that T^^(CN) 

~ 2. 8° K is obtained for several late-type stars which cannot 
maintain extensive H II regions. 

Although we lack exact knowledge of the physical conditions 
in most of the clouds observed, it is unlikely that electrons, 
probably the chief local contributor to the excitation in H I 
regions, are sufficiently abundant to affect T (CN) very much. 

One does not expect that the electron collision temperature 
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correction will exceed a few tenths of a degree, hence within 
the accuracy of the measurements one can probably equate 
T (12.64 mm) with T A (CN). 

This conclusion is susceptible to a direct test through 
observation of the CN X.2..6 mm radio lines: if the molecules are 
not in equilibrium with the background radiation then they will be 
observed in emission, by an amount dependent on the optical depth 
of the cloud. The molecular radiation temperature may be cal- 
culated if the cloud completely fills the antenna beam and has only 
a moderate optical depth, conditions quite probably fulfilled by 
the clouds considered here. If there is no background radiation, 
yet the excitation temperature is observed to be 2. 8° K, then the 
line radiation temperature of the strong transition F^F = 3/2 5/2 
-* \ 3/ 2 will be 


T d « 0. 25 t (° K) (23) 

(Thaddeus 1972) where is the optical depth of the line. For 
C °P h > where t q = 0. 34, T R 0. 08° K, barely detectable with 
existing receivers, but for BD+ 66° 1675, where t q ~ 5, this 
would be about 1. 2° K. Penzias, Jefferts, and Wilson ( 1972) 
have, in fact, made a search for the 3/2 5/2 - \ 3/2 line against 
BD+ 6 6° 1675 and find T R = 0. 036 ± 0. 037 °K, essentially a 

null result supporting the conclusion of the last paragraph. 
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It may be that this measurement, with the observed excitation 
temperature, offers the most convincing evidence of the validity 
of the CN thermometer. Two other stars where the CN optical 
depth is sufficient to make this sort of search worthwhile are 
13 Ceph and 20 Aql. 

While based on a small sample and hardly conclusive, the 
temperature invariance thus far observed fs highly indicative that 
the CN excitation results from a universal mechanism, and not 
from local interstellar conditions. 
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V. SUMMARY AND CONCLUSIONS 

A large number of high dispersion spectrograms of Q Oph 

and Q Per in the spectral interval AA3650 — 4350 A have been 

added together digitally to produce low noise, high resolution 

spectra. The improvement in signal-to-noise ratio 

has made possible the detection of six interstellar lines, of Ca I, 

13 + 

Fe I, K I and C H , in q Oph, their equivalent widths ranging 

o _ 

from 0. 3 to 2 mA, and three lines, of Ca I and Fe I lines of com- 
parable strength in Q Per. The syntheses have also provided new 
upper limits on as yet undetected lines relevant to atomic and 
molecular abundances — in particular to measurements 
of the intensity of the microwave background radiation at several 
millimeter and submillimeter wavelengths. 

o 

Ca I \4227 A, the most important of the newly detected lines, 
provides a direct determination of the ionization equilibrium and 

hence the electron density. In Q Oph this is found to be n = 0. 16 

-3 

— 0. 55 cm . In Q Per the existence of the Ca I line is 

marginal,, but if the suspected line is real one finds that n = 0. 044 

e 

-3 

— 0. 14 cm . The neutral hydrogen density is difficult to 
determine, but if the upper limit of Spitzer and Tomasko (1968) 

is adopted for the cosmic ray ionization rate of hydrogen, and the 
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depletion factor of carbon is assumed to be no greater than 10, 

-3 

then n(H I) in the Q Oph cloud lies between 50 and 700 cm ; 

-3 

in the Q Per cloud n(H I) lies between 10 and 80 cm 

In the Q Oph cloud sodium and potassium both appear at 
approximately their normal cosmic abundances relative to hydro- 
gen, but calcium, iron, and aluminum are underabundant by 
factors of 2500, 300, and > 100, respectively. The abundances 
relative to sodium in the Q Per cloud show substantially the same 
behavior. These results- support the idea that most of the iron, 
aluminum, and perhaps other elements are contained in the 
interstellar grains. 

o 

A new line at X4Z32 A in the Q Oph cloud is ascribed to 

C ^H + and yields a [C^/[C^] ratio of 81(+113, -16). This 

13 

is close to the terrestrial value 89, and suggests that C -rich 

stars have not contributed much to the interstellar medium (at 

least in the vicinity of the sun) and that the composition of the 

9 

interstellar medium has changed little in the past 5 x 10 years. 

The CN, CH, and CH^ interstellar bands have yielded a 

value of the intensity of the microwave background radiation at 

A.2. 64 mm and upper limits at 1. 32, 0. 559, and 0. 359 mm. The 

CN rotational temperature, T (\2. 64 mm), is found to be equal 

Jd 
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to 2. 75 ± 0. 10 °K for Q Oph and 2. 66 ± 0. 12 °K for Q Per. The 
short wavelength upper limits are consistent with a 2. 8° K 
blackbody spectrum and place constraints on the interpretation 
of recent rocket and balloon experiments. 

: Observations for the purpose of obtaining additional CN tem- 
peratures were made on two other early- type stars known to 
have CN absorption, and on two late B and five early A giants, 
two of which were found to contain CN. The rotational tempera- 
tures measured are all consistent with a 2. 8° K brightness 
temperature. 

It is interesting to note that no systematic factor of five 
increase in sensitivity in coude spectroscopy at the highest 
dispersion has been obtained since the advent of the coude 
spectrograph in the 1930's. Thus the results presented here 
encourage the belief that a number of new interstellar lines 
await discovery by the systematic application of new techniques. 

70 




64 


APPENDIX A 
DENSITOMETRY 

Wavelength and intensity calibration and spectral synthesis 
and filtering were done entirely digitally in the following way. In 
three successive densitometer scans of a given 10-inch plate the 
transmission T and displacment x at every 4 microns of plate 
travel were recorded digitally on magnetic tape for the two Fe 
comparison spectra and the stellar spectrum. In scanning the 
stellar spectrum the projected slit length of the densitometer 
was matched to that of the coude spectrograph so that all useful 
information was acquired in a single scan. Scans across the 
intensity calibration bars were next made at three different wave- 
lengths, and all data were then read into the large disk storage 
of the Goddard Institute IBM 360/95 computer system, with which 
all subsequent data analysis was done. 

By least squares fitting a third-order polynomial expansion 
of X as a function of x to about 25 Fe lines in both comparison 
spectra, an absolute wavelength calibration curve accurate to at 

o 

least 0. 004 A over the whole plate was obtained, and this 
was used to convert T to a function of X. The relative wave- 
length calibration from plate to plate was apparently even better 
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than 0. 004 A: on correcting for the earth's motion it was found 
that the wavelength assigned to a given strong interstellar line 
(e.g. H, K, CH + X4232, CH X4300) on different plates agreed to 
about 0. 002 A, which corresponds in radial velocity to only 0. 16 
km sec \ about 3% of the width of the sharpest interstellar lines 
in the -15 km sec ^ cloud. Thus no significant decrease in spectral 
resolution was to be expected from adding spectra together, and 
it was indeed found that those strong interstellar lines visible 
on a single plate showed no increase in width as successive 
spectra were folded into the final synthesis. • 


The intensity calibration function ’is that suggested by 
de Vaucouleurs (1968) for use in the threshold and toe regions 
of the characteristic curve: 



(Al) 


where I is the intensity, T and T^ are the spectrum and clear 

plate transmissions in arbitary units as measured by the densi- 

2 

tometer, and A and n are constants determined by a y fitting 
procedure. This function has the convenient property that when 
log [(T /T) - 1] is plotted against log I, a straight line results, 
simplifying the numerical calculations. Clearly, n may be 
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identified with the inverse of y, the slope of the linear part of 
this characteristic curve, so that this function will be appro- 
priate up to the point where the emulsion begins to saturate 
(D > 2 for Ila-O). The proportionality constant A has no actual 
importance since in the measurement of equivalent widths, only 
relative intensities of adjacent points in the spectrum are 
involved. In principle T could be measured, but in practice 
adjacency effects and the variation in fogging across some plates 
make this measurement often rather arbitrary. It was thought, 
then, that more precision could be obtained by introducing 
as a third parameter in the fit. The measurements of the cali- 
bration bars from the three scans across the plate were averaged 
together to obtain a single calibration function iather than being 
employed to determine a separate calibration for each region of 
the plate. This procedure is justified by the small variation of 
the low- density part of the Ila-O characteristic curve with wave- 
length in the region covered by the spectrograms (Latham 1969; 
Wright 1969). It is estimated that errors in equivalent width 
due to photometric reduction are under 4 percent. As a con- 
sistency check the strong H and K lines of Ca II were measured 
on single plates and the rms scatter of these equivalent widths 
was found to be 5 percent. This figure includes errors due both 
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to the intensity reduction and to placement of the intensity continuum. 
Finally, the results for strong lines from this synthesis agree quite 
well with those obtained by Herbig (see Table 1) with the same 
spectrograph. This comparison, of course, avoids the notorious 
discrepancy in equivalent widths measured at different observatories. 

Once the stellar spectrum is reduced to normalized intensity 
versus wavelength, it is corrected for the earth's orbital motion 
and converted by numerical interpolation to a standard form with 

o 

uniform spacing (0. 005 A). It then is added into a cumulative syn- 
thesis with a weight inversely proportional to the mean square 

o 

noise averaged over several 3 A intervals, following the usual ... 
prescription for maximizing the synthesis signal- to- noise ratio. 

Since the spectrograms were fairly uniform in quality, this refine- 
ment over equal weighting probably gained very little. 
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APPENDIX B 

STATISTICAL ANALYSIS OF UPPER LIMITS 

Since some of the more interesting results in this paper 
are upper limits on the strength of absorption lines, it becomes 
important to extract the maximum information possible from the 
synthesis. Here we present a simply statistical method based 
on the assumption that there is no probability for line emission, 
that improves to some extent the upper limits obtained from the 
usual two or three standard deviation criterion without loss of 
strictness. 

First consider the simple case of a rectangular absorption 
line one data point spacing in width which subtracts AI q from the 
continuum at one known wavelength and assume that the con- 
tinuum has a component of uncorrelated gaussian noise with 
standard deviation a. The object of the method is to find, with 
some arbitrary probability or strictness P, an upper limit s 

o 

to the true AI q » given o and the measured deviation at that 
wavelength AL The probability density is gaussian: 

p(Al|s) = (2tt a) exp [-(AI-s) 2 /2a 2 ] 

where s is the hypothetical true value for AI • For use in a 

o 
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hypothesis test, p(Al|s) should be multiplied by the probability 
density of s, which is assumed zero where s represents emission, 
but is unknown otherwise. Following the method first applied by 
Gauss we assume the probability density to be a constant where s 
would indicate absorption. The probability that s is contained in 
the interval [0, S] is then 

f* S ° 

1 p(Al|s)ds 
P (s )= tt 

J p(Al|s)ds 
o 


where the denominator provides the normalization. Inversion of 
this equation gives two cases: 


1) AI ^ 0 ("absorption") 


? _ 1 
s = AI + V2 a erf 
o 


2) AI > 0 ("emission") 


P - (l-P)erf 


(-ai/^7) 


s = AI + 
o 


erf” 1 


P + (1-P) erf 




As AI becomes large negative the expression reduces to the usual 
one, but as AI goes through zero to large positive, s^ approaches 
zero assymptotically. 


(B2) 


(B3a) 


(B3b) 
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The shape of the sharp interstellar lines is mainly determined 

by the instrumental profile of the spectrograph, which is to a good 

apprximation, gaussian. Since the plate grain noise is expected to 

be nearly white, the familiar theorem due to North that the optimal 

filter of a line in white noise is the line itself will be approximately 

correct. The filtering procedure adopted was to use the gaussian 

filter whose width maximized the signal-to-noise ratio of the 

+ ° 

moderately strong CH lines at \\3745 and 3958 A. The resultant 
filter width (full width at half maximum) was 0. 050 A compared to 
the 0. 058 A measured from the lines and led to about a factor 2 
improvement in signal-to-noise for these lines. 

In order to determine whether or not this approximate -v 

method was valid, we made a statistical analysis of the grain 
noise. The noise power spectrum of the synthesis was determined 
by Fourier analysis and found to be a superposition of two more 
or less gaussian distributions. The first, about five times wider 
than the transform of the sharpest interstellar lines, is the 
expected white grain noise already filtered to some extent by 
the densitometer analyzing slit. The second is of the same ampli- 
tude but has about the same width as the line transform: its 
origin is not completely understood. The general filter function 
appropriate to the problem at hand is the infinite memory filter' 
of Zadeh and Ragazzini (1952) which has the particularly simple 

, 77 

form: 
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CO 

qU) - J* 

- co 


A^ (uu) i(JU X _d«j 
$(ijd) 6 2tt 


where A3(uu) *- s the Fourier transform of the observed line profile, 
and i (uu) is the noise power as a function of spatial frequency. 
When $(uu) is a constant (white noise), q( X) reduces to the line 
profile. Thus the effect of the narrow noise distribution will be 
to broaden the filter transform and hence to narrow the filter, as 
was seen empirically above. Since the sum of the two noise dis- 
tributions is also very roughly gaussian, the filter will not differ 
much from the gaussian line shape except in width, and the func- 
tional distortion is estimated to be no more than the error 
involved in assuming the observed lines to be gaussian. Because 
of this conclusion, and because the major uncertainty in deter- 
mining upper limits is the placement of the continuum, which is 
to some extent discretionary, any refinement over the empirical 
best gaussian filter would have little effect in decreasing upper 
limits. 

The quantities AI and a in the filtered spectrum are 
determined in the following way. A least squares linear fit is 
made to a 1 to 2 A interval containing the expected wavelength 
of the line for which one desires an upper limit: since the line 

is not observed it cannot affect the continuum fit. The difference 
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AIj between the continuum and the measured intensity at the 
appropriate wavelength, and the rms error of the fit a are used 
in Equations (B3) to obtain an upper limit on line depth. The 
corresponding equivalent width upper limit is related to this by 
a constant of proportionally deduced from the weaker observed 
lines. The strictness criterion P = 0. 99 in Equations (B3) was 
chosen for a confidence level of 99 percent. 
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APPENDIX C 

STATISTICAL ERRORS IN EQUIVALENT WIDTH MEASUREMENT 


Once we have characterized the noise by its power spectrum, 

or the fourier inverse of the power spectrum , the autocorrelation 

function, it is not difficult to determine the standard deviation of the 

equivalent width integral. This will be the product of the range of 

integration and the rms error of the mean q„ontinuum computed over 

that interval. Let the integration range A comprise N points spaced 

2 

evenly apart, and a be the variance of the noise, normalized to 
unity continuum. For simplicity, the continuum is shifted to zero. 
The variance of the mean continuum is then given by 



N- 1 


2 

N 




R(k§X) 


(Cl) 


(Davenport and Root 1958, p. 80), where R(X) is the autocorrelation 

function of the noise. When the noise in uncorrelated over the sampling 

period, R(\) approaches a delta function and the second term makes 

2 

no contribution. The variance of the mean is then a /N, the result 

expected for the average of N statistically independent points. When 

the noise is highly correlated, however, R(\) is nearly a constant, 

2 2 

equal to a (since R(0) = a ), over the range of integration, and one 
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has a = a . Any single point measurement is as good as an 
average of adjacent highly correlated points. 


The noise autocorrelation function for the Q Oph synthesis is 

# 

to a good approximation 


R(X) = a 2 || 


exp 


-V /2 (0. 00'5 A) 


2 ]+ 7 exp[- 


X /Z(0. 029 A) 


(C2) 


o 

For N = 24, or A = 0. 120 A, more than 99 percent of the equivalent 
width is included in the integral. From equation (Cl) we have 
a M = 0. 40 a or = 48 a mi. 
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TABLE 1. INTERSTELLAR LINES IN Q OPHIUCHI IN 
THE REGION 

U 3650-4350 A 


Rest 

Wavelength 

o 

(A) 

Atom 

or 

Molecule 

f- value ^ 

This 

W* 

(mi) 

Work 

9 

V 

(km sec) 

Herbig 
W v 

o ” J- 

(mA ) (km sec ) 

3719. 935 

Fe I 

0. 041 h 

1. 9 

- 14. 9 " 

<3. 


■ 3745. 310 

CH* 


6. 0 

-15. 1 

7. 2 

-15. 1 

3779. 18 

CO + 

0. 0005(f 

<0. 20 


<3. 


3859. 913 

Fe I 

0. 023 J ' 

1. 0 

-15.0 

<3. 


3873. 369 a 

CN R(2) 


<0. 17 


< 3. 


3873. 998 

CN R( 1) 

0. 0069, 

2. 7 

-15. 0 

3.4 

-13. 8 

3874. 608 

CN R(0) 

0. 0173* 

8. 2 

-15. 2 

9.2 

-14. 8 

3875. 763 

CN P(l) 

0. 0058 

1.4 

-15.2 

^2. 

-14. 9 

3878. 768 

CH 

0. 0004/ 

2.0 

-14. 5 

3. 

-13.3 

3886.410 

CH 

0. 0014 1 

5. 7 

-15. 1 

5.9 

-14.4 

3888. 646 b 

He I 

_ _ , . 1X1 

0. 064 

<0. 25 


<3. 


3890. 213 

CH 

0. 00093 1 

3. 7 

-14. 7 

5.6 

-14. 3 

3933. 664 

Ca II 

0. 69 

34. 8 

-15. 1 

34-2 

-15- 0 




9-8 

-29. 1 

10.1 

-29. 0 

3944. 018 c 

A 1 I 

0. 115 

<0.26 


<3. 


3957. 700 

CH* 


14.8 

-15.7 

13. 3 

-15.4 

3758. 14 d 

C n H + 


<0. 25 


<3. 


3968.470 

Ca H 

0. 344 

20. 6 

-14. 8 

21. 3 

- 14. 2 




4. 7 

-28. 8 

5.4 

-28. 0 

3999. 08 

CO* 

0.00052^ 

<0. 23 


<3. 


4044. 145 b 

K I 

0.0061 

0. 8 

' -15. 3 11 

<3. 


4047. 213 

K I 

0. 0031 

0. 3 

-15.0 11 

<3. 


4119.48 

SiH 


<0. 22 


<3. 


4226. 728 

Ca I 

1. 75 

1. 3 

-15. 8 11 

<3. 


4229. 337® 

CH* R( 1) 


<0. 22 


<3. 


4232. 28 f 

c 13 h* 


0.4 

-15. 2 

<3. 


4232. 539 

CH* R(0) 


24. 3 

-15.0 

27.4 

-14. 0 

4250. 94 

CO + 

0. 00056 

<0. 28 


<3. 


4300. 321 

CH R 2 (l) 

0. 0052 1 

19.6 

-16.0 

20. 5 

- 14. 6 

4303. 947 g 

CH R x (l) 


<0. 28 


<3. 



"Except where noted, all rest wavelengths are from Herbig (1968). 

t Except where noted, all atomic f- values are as given by Wiese, Smith, and Miles (1969). 
^Upper limits represent a confidence level of 99 percent. 

J Heliocentric velocity of line centroid. 

HVelocity of line peak. 82 









Notes to Table 1. 


a ’ Jenkins and Wooldridge (1938). 

b ‘ Moore (1945). 

c. 

Burgess, Field, and Michie (I960). 

d. 

Bortolot and Thaddeus (1969). 

e. 

Douglas and Herzberg (1942). 

f. * 

' This work (see § lid). 

g " Gero (1941). 

b ' Bridges and Wiese (1970). 

1- Herbig (1948), Table 12. 

^‘Corliss and Bozman (1962), scaled to (h). 

^’Herbig (1968), Table 9. 

L jbid. , Table 8. 

m 


" Osterbrock (1964). 



TABLE 2. ATOMIC COLUMN DENSITIES IN THE 
C O p H - 15 KM SEC' 1 CLOUD 


atom 

uh 

N(cm ^ X 10 *) 

Fe I 

3719. 935 

3.8 


3859. 913 

3-3 

Ca I 

4226.728 

0.047 


* 


Ca II 

3933.664 




* 


3968.470 

4.8 

K I 

4044. 145 

9. 1 


4047. 213 

6.8 

A1 I 

3944.018 

<0.6 

He I 

388 . 646 

<0. 3 


-r 

Determined by doublet ratio method 
(Stromgren 1948), b= 1.92 km/sec. 









TABLE 3 


MOLECULAR COLUMN DENSITIES IN THE 


£ OPH - 15 KM SEC" 1 CLOUD 


Molecule. 

\(A) 

-2 -13 * 

N(cm x 10 ) 

CH - 

4300.321 

2.0 


3890. 213 

2.3 


3886.410 

2.4 


3878.768 

2.4 

CN 

3874.608 

0. 3 0 


3873. 998 

0.23 


3875.768 

0. 14 

SiH 

4119.48 

<0.04 

CO + 

- 4250.94 

<3.1 


3999.08 

<3.1 


3779. 18 

<3.2 


Calculated for b = 1. 36 km sec 

8S 








TABLE 4. ATOMIC PHOTOIONIZATION AND RECOMBINATION 


RATES IN THE - OPH - 15 KM SEC _1 CLOUD 


Atom 

/ ,„11, -1 
r (x 10 ) sec 

high field low field 



, in 11 * - 3 - 1 

a (x 10 ) cm sec 

Na I 

3. 96 

0.93 

<» 

0. 59 

A1 I 

81.4 

27.8 

0. 85 f 

K I 

10.4 

2.62 

0. 56* 

Ca I 

33.6 

9.6 

0. 60* 

Fe I 

44.2 

i. 9.6 

0. 60 § 


Seaton (1951) 

^Burgess, Field, and Michie (I960) 


assumed 
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TABLE 5. ATOMIC ABUNDANCES IN THE r CLOUD 
- 15 KM SEC" 1 CLOUD 










TABLE 6. INTERSTELLAR LINES IN Q PERSEI 
IN THE RANGE U3650 - 4350 


Rest 

Wavelength 

(A) 

Atom or Molecule 

W 

(ml) 

v* 

(km sec ) 

3719. 935 

Fe I 

1. 7 

12. 8 

3859.913 

Fe I 

1.7 

13.4 

3873. 369 

CN R(2) 

< 0. 31t 


3873. 998 

CN R( 1) 

3.2 

13. 0 

3874. 608 

CN R(0' 

11. 3 

12.8 

3875.763 

CN P(l) 

1. 3 

12.6 

3878.768 

CH 

1.6 

13. 5 

3886. 410 

CH 

4.9 

13.4 

3890. 213 

CH 

3.2 

13.7 



( 2. 2 

-5.8 

3933.664 

Ca II 

< 54. 5 

12.4 



l 5. 5 

20. 7 

3944.01 

A1 I 

< 0. 8t 


3957.700 

CH+ 

2.5 

12. 9 



( 2.6 

-5.3 

3968.470 

C a II 

<35. 5 

12. 8 



* 3.4 

21.0 

4226. 728 

Ca I 

0. 8 

12. 8 

4229. 337 

CH+ R( 1) 

< 0. 56f 


4232. 535 

CH+ R(0) 

2. 8 

12. 3 

4300. 321 

CH R 2 (l) 

15. 5 

12. 5 

4303. 947 

CH R 1) 

< 0.48 

J 


Heliocentric veLocity of line centroid. 


99 percent confidence level upper limits. 
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TABLE 8. MOLECULAR COLUMN DENSITIES IN 
THE + 13 KM SEC' 1 CLOUD IN r PER 


Molecule 

\ (A) 

N (cm -2 ) X 10' 13 * 

C H 

4300. 321 

1. 5 

' 

3890. 213 

2.0 


3886. 410 

2. 0 


3878. 768 

1.9 

C N 

3874.608 

0. 15 


3873. 998 

0. 10 


3875. 768 

0. 05 


Calculated for 


b + 1. 98 km sec 
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TABLE 9 


C 


Atom 


Na I 
A1 I 
Ca I 
Fe I 




rOMIC PHOTOIONIZATION RATES IN THE 
PER + 13 KM SEC' 1 CLOUD 



„ , ,- 11 . -1 
r (x 10 ) sec 

high field low field 



0.88 


73* 

23. 


27. 

8.5 


27 


9.3 








TABLE 10. TOTAL ATOMIC COLUMN DENSITIES IN THE 


c PER + 13 KM SEC' 1 CLOUD 


Atom 

N (cm 

Na 

in 15 

1.8 X 10 

Ca 

9.6 x 10 11 


t - 14 # 

Fe 

1.3 x 10 

- 

¥• 

„ 13 

A1 

<3. 2 x 10 


From a weighted mean of the two Fe lines.. 








TABLE 11. A COMPARISON OF THE ATOMIC ABUNDANCES RELATIVE TO 

i 

SODIUM OBSERVED AGAINST r OPH AND Q PER 




Fractional Abundance 

Atom 

e oph 

C Per 

Solar System 
(Cameron 1968) 

Na 

1.0 

1.0 

1.0 

Ca 

0.00079 

0. 00053 

1.2 

Fe 

0.070 

0.72 

14. 1 

A1 

<0.018 

<0.018 

1.3 
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TABLE 12. MOLECULAR EQUIVALENT WIDTHS AND OPTICAL 
DEPTHS IN THE Q OPH - 15 KM SEC* 1 CLOUD 


Molecule 

Line 

X . ( 1 ) 
air 

W(mA) t 

§ 

T 

CN 

R(0) 

3874.608 

8.27 ± 0. 10 

0. 338 ± 0. 003 


R(l) 

3873. 998 

2. 78 ± 0. 10 

0. 105 ± 0. 004 


P(l) 

3875. 763 

1.40 ± 0. 10 

0. 052 ± 0. 004 


R(2) 

3873. 369 

<0. 17 

<0.0062 

CH 

R 2 (D 

4300. 321 

19.6 

0. 719 


R X (D 

4303. 947 

<0. 28 

<0. 0081 

CH + 

R(0) 

4232. 535 

24. 1 

0.968 


R(l) 

4229. 337 

<0. 22 

<0. 0065 


Of the (0, 0) bands of the following electronic systems: 

CN B 2 Z - X 2 Z; CH A 2 A - X 2 Hj CH + A *11 - xhl . 

^Upper limits represent a confidence level of 99 percent; uncer- 
tainty on measurements is one standard deviation. 

§ ■■ 

Computed from the curve of growth for gaussian lines, using 
b = 1. 36 km sec . 
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TABLE 13. TEMPERATURE OF THE BACKGROUND RADIATION 
AT X2.64MM, FROM INTERSTELLAR CN IN Q OPH 

Line Ratio 


R(1)/R(0) P(1)/R(0) 


Tio(CN), uncorrected (°K) < 


2. 99 ± 0.07 2. 99 ±0.17 


Correction for curve of growth (°K) -0. 07 


-0.09 


Correction for electrons and H I (°K) 17 ± 0. 10 -0. 17 ±0. 10 


T (2. 64mm) (°K) 
B 


2.75. ± 0. 12 2.73 ± 0.20 


T n (2.64mm), weighted mean (°K) 
B 


2.75 ± 0. 10 


3 - 

Optical - depths derived from equivalent widths in Table 12, 
assuming, b = 2. 40 km sec * as measured from the unsmoothed 
spectrum in Figure 4. 

^Adopting b(CN) = 1. 36 km sec \ 

C Taking n^ = 0. 36 ± 0.19 cm ^ andn(HI) = 350± 350 cm 

^Relative weights 4: 1. The quoted uncertainty i3 one standard 
deviation. This includes no contribution for uncertainty in the 
curve of growth correction, which is felt to be at most com- 
parable to the other uncertainties. 
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TABLE 14. UPPER LIMITS ON BACKGROUND RADIATION FROM 
INTERSTELLAR MOLECULAR LINES IN Q OPH 


Molecule 

Line Ratio 

Wavelength 

(mm) 

Upper Limit 
on T b (°K) 

Upper Limit on 
IyXlO^^ferg cm -2 
sec - -^ster' * Hz' *) 

CN 

R(2)/R(l) 

1. 32 

3. 37 

0.72 

CH 

R 1 (l)/R 2 (l) 

0. 559 

5. 23 

1.66 

CH + 

R( 1)/R(0) 

0. 359 

7. 37 

3-75 
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TABLE 15. STARS OBSERVED 





















TABLE 16. CN. EQUIVALENT WIDTHS AND OPTICAL DEPTHS 


Star 

Line 

W(mA) v 

t 

T 

S Per 

R(0) 

11. 3 ±0.19 

0. 274 ± 0.005 


R(l) 

3. 2 ± 0. 19 

0.072 ± 0.005 

X Per 

R(0) 

14. 7 ± 1. 5 

P 

0.45 ±0.05 


R(l) 

■ 

• • 

• 

6.9 ± 1. 5 

0.19 ±0.05 

C Oph 

R(0) 

8. 27 ± 0. 10 

0. 338 ± 0. 003 


R(D 

1. 40 ± O.'lO 

0. 105 ± 0. 004 

20 Aql 

R(0) 

9. 3 ± 0. 9 

0. 27 ± 0.03 


R(D 

3. 5 ± 0. 9 

0.095 ± 0.025 

HR7573 

- R(0) 

5.0 

0.14 


R(D 

<1.9 

<0.05 

13 Ceph 

R(0) 

19.3 ±1.7 

0.62 ± 0.07 


R(l) 

7.4 ±1.7 

0.21 ±0.05 


^Uncertainty on measurements is one standard deviation; upper 
limit represents a confidence level of 99 percent. 


^Computed with b = 1.98 km sec ^ for CPer, 1. 36 km sec ^ for 
Q Oph, and 1.67 km sec ^ for all other stars. 

98 



















REFERENCES 


Adams, W. S. 1941, Ap. J. , 93, 11. 

, 1949, Ap. J. , 109, 354. 

Allison, A. C. , and Dalgarno, A. 1971, Astron. and Astrophys. , 

13, 331. 

Alpher, R. A. , Bethe, H. , and Gamow, G. 1948, Phys. Rev. , / 

73, 803. 

Alpher, R. A. , and Herman, R. C. 1949, Phys. Rev. 75 , 1089. 

, 1951, Phys. Rev. , 84, 60. 

Augason, G. C. , and Herbig, G. H. 1967, Ap. J. , 150 , 729. 

Beals, C. S. , and Oke, J. B. 1953, M. N. , 11 3, 530. 

Bennett, R. G. , and Dalby, F. W. 1962, J. Chem Phys. , 36 , 399. 

Blair, A. G. , Beery, J. G. , Edeskuty, F. , Hiebert, R. D. , 

Shipley, J. P. , and Williamson, K. D. , Jr. 1971, Phys. 

Rev. Letters, 27, 1154. 

Bless, R. C. , and Savage, B. D. 1970, in Ultraviolet Stellar Spectra 

and Related Ground Based Observations (I. A. U. Symp. No. 36), 
ed. L. Houziaux and H. E. Butler (Dordrecht: Reidel) 

Borgman, J. , and Blaauw, A. 1963, B. A. N. , 17, 358. 

Bortolot, V. J. , Jr., Clauser, J. F. , and Thaddeus, P. 1969, 

Phys. Rev. Letters, 22, 307. 

Bortolot, V. J. , Jr. , and Thaddeus, P. 1969, Ap. J. (Letters), 155, LI 7. 

Burgess, A. , Field, G. B. , and Michie, R. W. I960, Ap. J. , 131, 529. 

100 



Cameron, A. G. W. 1968, in Origin and Distribution of the Elements , 
ed. L. H. Ahrens (Oxford: Pergamon). 

Carruthers, G. R. 1970, Ap. J. (Letters), 16 1 , L81. 

Clauser, J. F. 1970, unpublished Ph. D. dissertation, Dept, of 
Physics, Columbia. 

Clauser, J. F., and Thaddeus, P. 1972, Proc, 3rd International 
Conf. on Relativistic Astrophysics , (to be published). 

Corliss, C. H. , and Bozman, W. R. 1962, Experimental Transition 
Probabilities for Spectral Lines of Seventy Elements; 

NBS Monograph 53. 

Davenport, W. B. , and Root, W. L. 1958, An Introduction to the 

Theory of Random Signals and Noise (New York: McGraw Hill), 
de Vaucouleurs, G. 1968, Applied Optics, 7 _, 1513. 

Dicke, R. H. , Peebles, P. J. E. , Roll, P. G. , and Wilkinson, D. T. 
1965, Ap. J. , _U2, 414. 

Ditchburn, R. W. , and. Hudson, R. D. I960, Proc. Roy. Soc. 

London (A), 256 , 53. 

Douglas, A. E. , and Herzberg, G. 1942, Can. J. Res. , 20A , 71. 
Field, G. B. , Goldsmith, D. W. , and Habing, H. J. 1969, Ap. J. 
(Letters), 155 , L149. 

Field, G. B. , and Hitchcock, J. L. 1966, Phys. Rev. Letters, 16, 817 
Fink, E. H. , and Welge, K. H. 1967, J. Chem. Phys. , 46, 4315. 
Fowler, W. A., Caughlan, G. R. , and Zimmerman, B. A. 

1967, Ann. Rev. Astr. and Astrophys. , jj, 525. 


101 





Frisch, P. 1972, Ap. J. , P73, 301. 


Gamow, G. 1948, Phys. Rev. , _74, 505. 

Gero, L. 1941, Z. Physik, _U8, 27. 

Goldsmith, D. W. , Habing, H. J. , and Field, G. B. 1969, Ap. J. , 
158 , 173. 

Habing, H. J. 1968, B. A. N. , 19, 421. 

. 1969, B. A. N. , 20, I 77 . 

Hegyi, D. , Traub, W. , and Carleton, N. 1972 (to be published). 
Helstrom, C. W. 1968, Statistical Theory <sf Signal Detection 
(New York: Pergamon). 

Herbig, G. H. 1968, Z. fur Astrophysik, ^_8, 243. 

Herzberg, G. 1959, The Spectra of Diatomic Molecules 
(New York: McGraw-Hill). 

Hickock, F. R. , and Morton, D. C. 1968, Ap. J. , 152 , 203. 
Hjellming, R. M. , Gordon, C. P. , and Gordon, K. J. 1969, 
Astron. and Astrophys. , 2, 202. 

Hobbs, L. M. 1965, Ap. J. , 142, 160. 

. 1969, Ap. J. , 157, 135. 

Hobbs, L. M. , and Zuckerman, B. 1972, Ap. J. , 17 1 , 17. 

Houck, J. R. , and Harwit, M. O. 1969, Ap. J. (Letters), 

157, L45. 

Howard, W. E. , Wentzel, D. G. , and McGee, R. X. 1963, 


Ap. J. , 138, 988. 


102 



Hudson, R. D. , and Carter, V. L. 1965, Phys. Rev. , 139 , A1426. 


. 1967a, Ap. J. , 149, 229. 

. 1967b, T. Opt. Soc. Am. , _57, 651. 

Jefferts, K. B. , Penzias, A. A., and Wilson, R. W. 1970, 

Ap. J. (Letters), 16 1 , L87. 

Jenkins, F. A. , and Wooldridge, D. E. 1938, Phys. Rev. , _53, 137. 
Ap. J. (Letters), 16 1 , L87. 

Kelly, H. , and Ron, A. 1971, Phys. Rev. 'Letters, Zb_, 1359. 

Kozlov, M. G. , Nikonova, E. I. , and Startsev, G. P. 1966, 

Optics and Spectr. , 2_1, 298 (USSR). 

Latham, D. W. 1969, Bull. A. A. S. , _1, 147. 

Lillie, C. F. 1968, unpublished Ph. D. dissertation, 

Dept, of Astronomy, Univ. of Wise. 

Lynds, B. 1969, P. A. S. P. , 81, 496. 

Mather, J. C. , Werner, M. W. , and Richards, P. L. 1971, 

Ap. J. (Letters), 170 , L59. 

McKellar, A. 1951, Publ. Dominion Astrophys. Observatory, 
Victoria, B. C. , _7, 251. 

. 1965, in Stellar Atmospheres, ed. J. L. Greenstein, 

(Chicago: Univ. of Chicago). 

Merrill, P. W. , Sanford, R. F. , Wilson, O. C. , and Burwell, C. G. 
1937, Ap. J. , 86, 274. 


103 



Moore, C. E. 1945, A Multiplet Table of Astrophysical Interest 


(Princeton: Princeton Observatory). 

Muehlner, D. , and Weiss, R. 1970, Phys. Rev. Letters, _24, 742. 

. 1972 (£o be published). 

Munch, G. 1957, Ap. J. , 125, 42. 

.1964, ibid., 140, 107. 

Newsom, G. H. 1966, Proc. Phys. Soc. , _87, 975. 

Osterbrock, D. E. 1964, Ann. Rev. Astr. and Astrophys. , _2, 95. 
Peimbert, M. 1968, Bull. Obs. Tonantzintla y Tacybaya, 4, 233. 
Penzias, A. A., Jefferts, K. B. , and Wilson, R. W. 1972, 

Phys. Rev. Letters, 28 , 772. 

Penzias, A. A. , and Wilson, R. W. 1965, Ap. J. , 142 , 419. 

Pipher, J. L. , Houck, J. R. , Jones, B. W. , and Harwit, M. O. 
1971, Nature, 2LL, 375. 

Rees, M. J. , Sciama, D. W. , and Stobbs, S. H. 1968, Ap. Letters, 
2, 243. 

Rogers, A. E. E. , and Barrett, A. H. 1966, A. J. , 71 , 868. 
Rogerson, J. B. , Spitzer, L. , Jr. , and Bahng, J. D. 1959, 

Ap. J. , 130, 991. 

Seaton, M. J. 1951, M. N. , 111 , 368. 

1959, ibid., 119, 81. 

Scherb, F. 1968, Ap. J. (Letters), 153 , L55. 

104 



Shane, C. D. , and Wirtanen, C. A. 1967, Publ. Lick Obs. , 

22 , Pt. 1. 

Shivanandan, K. , Houck, J. R. , and Harwit, M. O. 1968, 

Phys. Rev. Letters, 21 , 1460. 

Silk, J. , and Werner, M. W. 1969, Ap. J. , 158 , 185. 

Smith, A. M. 1972 (to be published). 

Smith, A. M. , and Stecher, T. P. 1971, Ap. J. (Letters), 164 , L43. 
Spitzer, L., Jr. 1948, Ap. J. , 108 , 276. „ 

. 1968, Diffuse Matter in Space (New York: Wiley). 

Spitzer, L. , Jr. and Tomasko, M. G. 1968, Ap. J. , 152 , 971. 
Stecher, T. P. 1968, Bull. Astr. Inst. Czech. , JL9, 279. 

Stromgren, B. 1948, Ap. J. , 108 , 242. 

Thaddeus, P. 1972, Ann. Rev. Astr. and Astrophys., (to be published). 
Thaddeus, P. , and Clauser, J. F. 1966, Phys. Rev. Letters, 

16, 8 19. ~ ~ 

Van Citters, G. W. , and Morton, D. C. 1970, Ap. J. , 16 1 , 695. 
Whitford, A. E. 1958, A. J. , _63, 201. 

White, R. E. 1972 (to be published). 

Wiese, W. L. , Smith, M. W. , and Miles, B. M. 1969, Atomic 

Transiti on Probabilities, Volume H: Sodium through Calcium, 
NSRDS-NBS 22. 

Wilson, O. C. 1948, P. A. S. P. , 60, 198. 

Witt, A. N. , and Lillie, C. F. 197.1, (to be published). 

105 



Wright, K. O. 1969, Bull. A. A. S. , _1. 14 2. 

Zadeh, L. A. , and Ragazzini, J. R. 1952, Proc. I. R. E 
40. 1223. 

Zimmermann, H. 1964, Astron. Nachr. , 288, 99. 




106 






Figure 



VELOCITY (KM/SEC) 


2. The atomic line profiles observed against Q Oph in the 
interval \\3650 — 4350 A plotted against heliocentric 
velocity. 108 
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Figure 3 


The molecular line profiles observed against Q Oph in the 
interval \\3650 — 4350 A plotted against heliocentric 
velocity. 
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Figure 3. (cont'd) 
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Figure 4. The \387 5 A CN band in Q Oph. Trace (a) is the unfilter- 

o 

ed synthesis, (b) is filtered with a 0. 050 A FWHM gaussian 
and (c) is a magnification of (b). The vertical line marks 
the position of R(2) and the short line through the magni- 
fied spectrum is the best straight line continuum fit to 
the interval it subtends. Ill 










Figure 7. 


Figure 8. 


The high and low radiation fields adopted at the 
position of the Q Oph-15 km sec ^ cloud. The high 
field is that of Q Oph at a distance of 15 pc plus the 
galactic field of Habing (1968), both attenuated by the 
cloud itself. The low field is the attenuated galactic 
field alone. 

The adopted photoionization cross-section of Ca I. 

o 

From \A2028 to 1660 A, the recommendations of 
Hudson and Carter (1967a) have been followed. The 
cross-sections of Newsom (1966) have been adopted 
after multiplication by 1. 5, except at the peaks of 
the strong AA1886 and 1765 A autoionizing lines, 
where Hudson and Carter's own values were accepted. 
In the range W1660 to 1080 A, 1. 5 times the cross- 
sections of Ditchburn and Huson (I960) have been 

o 

used, while below \1080 A the constant value 
a = 1. 0 Mb has been assumed. 

A.' 


114 



ABSORPTION CROSS SECTION Mb 





ELECTRON DENSITY (CM - 



Figure 9. The electron density as a function of neutral hydrogen 

density and trace element abundance. The horizontal 

- 3 

lines at n^ = 0. 16 and 0. 55 cm represent the limits 

- 1 

derived for the Q Oph - 15 km sec cloud. 
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Figure 11. The molecular line profiles observed against 

in the interval \\3650 — 4350 A, plotted against helio 


centric velocity. 
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c r** 




P(l) 



3875 3876 


band in Q Per. 





Figure 14. Direct and indirect measurements of the microwave back- 
ground radiation. Open circles are the direct ground- 
based measurements (see Thaddeus 1972 for references), 
closed circles are the indirect measurements obtained 
from the CN, CH and CH + interstellar bands in Q Oph 
(Tables 13 and 14), C represents the most recent Cornell 
rocket results, and M the most recent MIT balloon results 
after correction for atmospheric emission. 












